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Preface 


Currently,  research  is  being  conducted  by  the  Air  Force 
Wright  Aeronautical  Laboratories  to  replace  copper  wire  data 
links  on  board  aircraft  with  an  optical  fiber  system.  The 
optical  fibers  are  to  be  used  to  transmit  information  from 
sensors,  which  are  placed  throughout  the  aircraft,  to  a  cen¬ 
tral  receiver.  Because  of  the  attenuation  in  the  cable  and 
the  optical  to  electrical  conversion  process,  the  signal  in¬ 
curs  variations  causing  the  detection  accuracy  to  decrease. 
This  report  is  to  analyze  the  effects  of  adding  an  automatic 
gain  control  (AGC)  circuit  to  provide  amplification  of  the 
signal  and  control  to  the  signal  variations. 

I  would  like  to  thank  my  advisor,  LtCol  R.  J.  Carpinella, 
of  the  Air  Force  Institute  of  Technology,  for  his  guidance 
towards  the  completion  of  this  study.  I  would  also  like  to 
express  my  gratitude  to  my  fellow  classmates,  family,  and  my 
fiancee,  Lynette  Pray,  for  their  support  and  inspiration 
throughout  the  extent  of  this  study.  And  finally,  I  wish  to 
acknowledge  my  thanks  to  my  typist,  Mrs  Ruth  Walker. 
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This  paper  is  to  analyze  the  effects  of  an  automatic  gain 
control  (AGC)  in  an  optical  receiver  system.  The  intent  of 
using  an  AGC  is  to  improve  the  performance  of  the  receiver. 

The  received  signal  is  Manchester  encoded  and  contains 
fluctuations  caused  by  transmission  attenuation.  In  the  con¬ 
version  from  an  optical  signal  to  an  electrical  signal,  the 
photodetector  produces  shot  noise.  The  purpose  of  the  AGC 
will  be  to  amplify  the  photodetector  signal,  control  signal 
fluctuations,  reduce  shot  noise,  and  improve  detection  ac¬ 
curacy  . 

Using  this  criteria,  the  effects  of  the  AGC  in  the 
receiver  was  examined.  But  because  of  the  non-linearities 
associated  with  an  AGC,  an  analytical  analysis  can  not  be  con¬ 
ducted.  Therefore,  the  receiver  was  simulated  on  the  com¬ 
puter.  This  paper  describes  the  models  and  algorithms  used  to 
analyze  the  receiver.  Through  the  simulations,  the  receiver’s 
performance  with  an  AGC  is  compared  to  the  performance  without 
an  AGC. 

Because  of  the  limited  time  for  the  analysis,  only  three 

specific  AGC's  were  examined.  The  first  AGC  studied  is  clas- 

< _ 

/ 

sified  as  an  instantaneous  AGC  (IAGC)  followed  by  an  analysis 
of  a  common  AGC.  The  final  AGC  to  be  investigated  is  a  dual 
AGC  which  is  an  AGC  cascaded  to  an  IAGC. 
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IMPROVED  SIGNAL  CONTROL:  AN  ANALYSIS  OF  THE  EFFECTS  OF 
AUTOMATIC  GAIN  CONTROL  FOR  OPTICAL  SIGNAL  DETECTION 


I.  Introduction 


Background 

In  recent  years  the  application  of  optical  fibers  in  com¬ 
munication  systems  has  increased  significantly.  Technological 
developments  have  enabled  low  loss  optical  fibers  to  replace 
copper  wire  communication  links.  Improvements  were  made  in 
laser  diodes  and  LED’s  to  achieve  higher  optical  power  out¬ 
puts.  Because  of  the  higher  output  powers  and  low  loss 
cables,  the  distance  between  repeaters  can  be  increased, 
thereby,  reducing  construction  cost.  Construction  costs  were 
also  decreased  and  the  construction  simplified  by  improved 
techniques  and  hardware  for  the  interfacing  and  the  splicing 
of  optical  fibers.  These  are  just  a  few  of  the  advancements 
which  have  contributed  to  the  feasibility  of  an  optical  fiber 
communication  system  and  to  its  resultant  growth. 

The  practicality  of  an  optical  fiber  communication  system 
explains  the  current  growth  and  increased  usage,  but  the  in¬ 
herent  advantages  of  the  system  prompted  the  initial  interest 
and  research.  The  most  important  attribute  is  that  more  in¬ 
formation  can  be  transmitted  through  optical  fibers  than 
through  copper  wire.  Currently,  Gigabit  data  rates  have  been 
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obtained  by  using  single  mode  fibers.  Optical  cables  are  also 
immune  to  electromagnetic  radiation.  No  special  shielding  is 
required  to  protect  them  against  high  electric  or  magnetic 
fields.  The  fibers  are  lightweight  and  pose  no  hazard  to  fire 
or  electrical  shock.  Because  of  these  advantages  and  the 
technological  developments,  optical  fibers  have  become  a  cost- 
effective  replacement  for  copper  wire;  and  optical  fiber  com¬ 
munication  systems  are  being  implemented  in  many  new  projects. 

Problem 

One  such  project  is  to  implement  an  optical  fiber  data 
link  on  board  an  aircraft.  The  optical  cables  would  carry 
information  from  several  sensors  to  a  central  receiver.  The 
sensors  are  placed  at  various  positions  throughout  the  air¬ 
craft.  Each  sensor  will  transmit  a  Manchester  encoded  optical 
signal  via  an  optical  fiber.  Because  the  signals  must  travel 
different  path  lengths  to  reach  the  receiver,  the  attenuation 
losses  will  cause  the  optical  power  level  to  fluctuate  between 
the  signals.  At  the  receiver  several  bits  of  information  from 
each  signal  will  be  observed  for  one  frame  period  (Tf).  Thus, 
the  received  signal  (see  Figure  l)  will  be  a  Manchester  en¬ 
coded  optical  signal  with  fluctuations  in  the  power  level 
depending  on  which  sensor  is  being  observed. 
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Figure  1.  Received  Optical  Signal 


The  received  optical  signal  is  then  converted  into  an 
electrical  signal  by  a  photodetector-  Because  of  the  fluctua¬ 
tions  in  the  received  signal,  the  detector’s  output  signal 
will  vary  over  a  wide  dynamic  range.  The  photodetector  signal 
will  also  incur  variations  caused  by  the  dark  current  and  shot 
noise  which  are  inherent  to  the  photodetectors.  The  shot 
noise  is  used  to  describe  the  randomness  associated  with  the 
detector's  generation  of  an  electrical  signal.  The  dark  cur¬ 
rent  is  a  product  of  thermally  emitted  electrons  flowing  from 
the  detector.  The  dark  current,  shot  noise  and  signal  fluc¬ 
tuations  complicate  the  process  of  reconstructing  the  informa¬ 
tion  sent. 

For  the  particular  signal  being  studied,  the  information 
is  digitally  encoded  and  is  extracted  by  a  decision  circuit. 
The  decision  circuitry  determines  whether  a  binary  "one”  or 
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"zero"  vas  sent  based  on  ‘•he  photodetector's  output.  To 
increase  the  accuracy  of  the  decision  process,  an  automatic 
gain  control  (AGC)  circuit  will  be  placed  in  between  the 
photodetector  and  the  decision  circuitry.  The  purpose  of  the 
AGC  will  be  to  amplify  the  detector  signal  and  to  auto¬ 
matically  adjust  its  gain  to  reduce  signal  variations, 
thereby,  hopefully  improving  the  accuracy  of  the  reconstruc¬ 
tion. 

In  this  paper  a  detailed  analysis  will  be  conducted  on 
the  receiver  to  investigate  the  effects  of  an  AGC  in  an  op¬ 
tical  receiver.  Because  of  the  non-linearities  associated 
with  an  AGC,  an  analytical  analysis  of  the  receiver  becomes 
very  complicated.  Therefore,  an  empirical  analysis  will  be 
conducted  using  a  computer  to  simulate  the  receiver.  Through 
the  use  of  the  computer,  the  effects  of  the  AGC  on  the  detec¬ 
tor  signal,  the  signal  to  noise  ratio,  and  the  probability  of 
bit  errors  will  be  examined. 

Scope 

Because  of  the  limited  time  to  conduct  an  analysis,  only 
three  AGC’s  will  be  studied.  The  first  AGC  to  be  analyzed 
will  be  an  instantaneous  AGC  (IAGC),  followed  by  an  analysis 
of  a  common  AGC.  Both  AGC's  employ  a  single  loop  feedback 
configuration  to  provide  the  gain  control.  They  both  also 
have  a  gain  characteristic  which  is  an  exponential  function  of 
the  control  voltage.  The  final  AGC  to  be  studied  will  be  a 


dual  AGC,  which  is  a  cascaded  combination  of  an  IAGC  and  an 


AGC 


The  photodetector  will  be  a  non-amplifying  photodiode, 
such  as  a  p-n  or  a  p-i-n  diode.  The  effects  of  avalanche 
photodiodes  will  not  be  considered  in  this  analysis. 


Assumptions 

The  photodetector  will  be  illuminated  solely  by  the  Man¬ 
chester  encoded  optical  signal.  The  detector  will  be  com¬ 
pletely  enclosed  to  eliminate  the  effects  from  the  background 
radiation.  The  received  optical  signal  will  be  undi3torted 
except  for  the  attenuation  caused  by  the  transmission  through 
the  optical  fiber.  The  signal  is  also  assumed  to  be  spatially 


coherent,  so  that  the  optical  power  can  be  written  as 


P(t)  -  ffk  U(r,t)  dr  -  p(t)  A  (1) 

where 

U(r,t)  ■  spatial  and  time  representation  of  the  signal 
power  density 

p(t)  *  time  varying  power  density 
A  *  area  of  the  detector 

All  noises  will  be  neglected  except  for  the  shot  noise  and 
dark  current.  The  timing  of  the  electrical  circuit  will  be 
exactly  synchronized  with  the  optical  signal. 


Approach  and  Presents  tion 

Before  any  analysis  of  the  receiver  is  conducted,  the 
performance  parameters  of  the  receiver  will  be  defined.  Based 
on  the  photodetector  output,  an  optimum  decision  rule  is 
developed  to  minimize  the  probability  of  a  bit  error.  The 
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effects  that  the  AGC  has  on  the  probability  on  error  and  the 
signal  fluctuations  will  be  used  to  assess  the  AGC’s 
performance.  The  objective  of  adding  an  AGC  will  be  to  am¬ 
plify  the  detector  output,  eliminate  the  signal  fluctuations 
and  to  decrease  the  probability  of  a  bit  error. 

With  these  objectives  in  mind,  the  mathematical  models 
are  developed  to  describe  the  receiver.  These  models  need  to 
represent  the  physical  devices  as  accurately  as  possible  to 
obtain  valid  results.  From  the  mathematical  models  computer 
algorithms  will  be  developed.  Some  simple  simulations  of  the 
photodetector  and  AGC  will  be  conducted  to  test  the  validity 
of  the  computer  programs. 

Using  these  computer  programs,  the  effects  of  the  three 
specific  AGC's  on  the  photodetector  signal  will  be  analyzed. 
The  AGC's  response  to  the  detector's  output  will  be  obtained 
by  plotting  the  first  and  second  order  statistics  of  the  AGCs’ 
output.  The  effects  of  the  AGC  on  the  probability  of  error 
will  also  be  examined.  From  this  analysis  conclusions  can  be 
drawn  as  to  the  effectiveness  of  an  AGC  in  an  optical 
receiver . 
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IT.  Performance  Parameters 


The  purpose  of  a  receiver  is  to  reconstruct  as  accurately 
as  possible  the  information  which  was  transmitted.  Unfor¬ 
tunately,  the  signal  becomes  distorted  during  the  transmission 
and  the  reception.  The  distortion  introduces  error  in  the 
reconstruction  of  the  information  sent.  For  the  digitally 
encoded  signals,  the  information  is  transmitted  in  binary 
form.  Each  single  piece  of  binary  information  sent  is  re¬ 
ferred  to  as  a  bit  and  is  represented  by  a  "one"  or  "zero". 
Distortion  of  a  digital  signal  causes  incorrect  decisions  or 
bit  errors  to  be  made.  In  the  following  section  the 
probability  of  a  bit  error  is  minimized  through  the  develop¬ 
ment  of  an  optimal  decision  rule. 

In  the  last  section  of  this  chapter  the  objectives  of 
adding  an  AGC  will  be  discussed.  These  objectives  will 
specify  the  parameters  to  assess  the  performance  of  the  AGC's. 

Optimum  Decision  Criteria 

For  the  digital  signals  which  are  being  analyzed,  the 
binary  information  is  determined  by  the  position  of  the  pulse. 
This  type  of  digital  encryption  is  called  pulse  position 
modulation  (PPM).  A  binary  "one"  will  be  assigned  to  the  op¬ 
tical  pulses  which  occur  during  the  first  half  of  a  bit  period 
(Tg) ,  and  a  binary  "zero"  to  the  pulses  which  occur  in  the 
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second  half  of  the  bit  period. 

At  the  receiver  these  optical  pulses  are  converted  into 
an  electrical  signal  by  a  photodetector.  The  electrical  sig¬ 
nal,  that  is  generated  by  the  detector,  is  no  longer  a  deter¬ 
ministic  signal;  but  becomes  a  weighted  Poisson  counting 
process  (will  be  proven  in  the  next  chapter).  Because  the 
signal  is  now  a  random  process,  errors  will  occur  in  the  deci¬ 
sion  process.  Fortunately,  the  rate  of  the  counting  process 
is  proportional  to  the  power  density  of  the  optical  signal. 
Thus,  the  binary  signals  can  be  represented  by  the  rates  of 
the  counting  processes  (see  Figure  2).  An  optimum  decision 
rule  will  be  developed  to  extract  this  deterministic  signal 
from  the  Poisson  counting  process,  N(t). 


Xg  -  rate  of  photon  to 
.  electron  conversions 
Ap  -  rate  of  thermally 
emitted  electrons 


Figure  2.  Deterministic  Binary  Signal  of  the  Photodetector 


To  obtain  a  decision  rule,  the  Poisson  counting  process 


is  to  be  equivalently  written  in  a  vector  form.  This  vec- 
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torization  is  accomplished  by  sampling  the  counting  process, 
N(t),  at  a  rate  defined  by: 


where 

*  the  bandwidth  of  the  counting  process 
This  sampling  scheme  produces  a  count  vector  given  by: 


(2) 


k  =  (kj  ,k2 . k j)  (3) 

where 

kj  "  *(tj)  "  5  for  ,2 . J 

J  -  Tf 

s 

The  probability  of  a  count  k^  occurring  is  defined  by  (Ref 
3:205+) : 

k  .  X  .  .At 

Pr(kj/X1j)  =  t(XijAt)  i  e  13  3  /  kj!  (4) 

where 


At 

X 


t . 
J 


Vi 


th 


ij 


rate  of  the  i  signal  during  the  interval 

VVi 


Since  the  probabilities  of  each  sample  count  are  independent, 
the  probability  of  the  count  vector  can  be  written  as: 


_  J 

Pr(k/X  )  -  It  [(X 

j-1 

Based  on  these  probabilities  a  decision  rule  can  be  developed. 

The  optimum  decision  rule  for  minimum  probability  of  er¬ 
ror  can  be  derived  from  the  maximum  a  posteriori  (MAP)  deci- 


U 


At)  j 


X.  .At 
.  1J  ]  / 


k  ’ 

y 


(5) 


9 


3  i  on  criteria. 


The  MAP  decision  criteria  is  defined  as: 


1 

A(r)  -  Tp(X  /k)/P(x  /k) ]  <  1  (6) 

1  0  D 

0 

where 

A(r)  =  likelihood  ratio 

P(X^/k)  ■  probability  that  the  rate  is  X^  given  the 
observed  vector  k 

By  applying  Bayes  rule,  the  MAP  criteria  can  be  equivalently 
written  as: 

f Pr(k/X  )  Pr(X  ) ]  /  Pr(k)  >1 

A  (  r )  =  - 1 - —  <  1  (7) 

[ Pr(k/X  Q)  Pr(X0)]  /  pr(k)  D 

where  ° 

Prfk/X^)  -  defined  by  eqn  (5) 

Pr(X^)  *  priori  probabilities  of  the  given  rate  X^ 

Pr(k)  *  probability  of  the  vector  k  occurring 
Assuming  equal  priori  probabilities  and  substituting  eqn  (5) 

into  eqn  (7)  produces: 


A  (r) 


J  k 

n  T(x  At) 
j-1  Oj 

Since  the  logarithm  is 
equivalent  MAP  rule  is: 


A 


-X  At 

•  1  '  V 

-A. .At 

e  ^  1  /  V 


a  monotonic  increasing  function, 


(8) 


an 


1  0 


1QA(  r) 


£  [kj  l«iu4t  -  - 

J  ^ 

J 

E  [k .  In  A  At 

j-t  3  01 


>Q  4t]  <'  0 
°0 


(9) 


By  separating  the  summations  and  correctly  defining  the  limits 
produces  constant  rates  for  each  of  the  smaller  summations. 
This  causes  eqn  (9)  to  be  rewritten  as: 


J/2  J 

E  [k.  Xon At  -  XonAt]  +  E  [k  In  X.At  -  X  At] 

j  =  1  3  SD  SD  j  = J/2  +  1  3  D  D 

J/2  J  ®1 

E  [k.  X-r.  At  +  Xfi  At  ]  -  E  [k.  lnX  nAt  -  X<,_At]  <  0  (lO) 

j-1  3  D  D  J-J/2M  0  SD  SD  Dq 

By  rearranging  the  terms  and  simplifying,  eqn  (10)  reduces  to: 


J/2  J  D. 

2k-rk.ro  (11) 

j  =  1  3  j  = J/2  + 1  J  Dq 

Thus,  the  optimum  decision  rule  compares  the  count  at  the  end 
of  the  first  half  of  a  bit  period  with  the  counts  which  occur¬ 
red  in  the  second  half  of  the  bit  period.  Writing  this  in 
terms  of  the  actual  counting  process  produces: 

D1 

N(T/2)  %  N(T)  -  N(T/2)  (12) 

D0 

The  probability  of  error  is  given  by: 


Pr (  error)  -  Pr ( N(  T/2)  >N(  T) -N(  T/2  )  /  Xq)  P(  Xq )  + 
1/2  Pr(»(T/2)-N(T)-N(T/2)/ Xq)P(  Xq)  ♦ 
Pr(N(T)-N(T/2)>N(T/2)/  X1  )P(  X1  )  ♦ 
1/2  Pr(N(T)-H(T/2)-N(T/2)/ X1  )P(  X1  ) 

1  1 


(13) 


Substituting  in  the  Poisson  probability  densities  eqn  (13) 
becomes  (Ref  3:220): 


“  i.  >  -XnT/2 

Pr(  error)  =  E  [  1  -  I  [x^T/P)*  e  ]  /  k: !  ] 

i=0  k=0 

,  -X„T/2 

[(XSDT/2)X  e  S  ] / i !  ]  (14) 

According  to  the  MAP  decision  criteria,  the  decision  rule 
produces  the  minimum  probability  of  error  (Ref  6:42).  The 
decision  rule  was  developed  assuming  constant  optical  power 
densities  of  all  the  received  pulses.  But  the  decision  rule 
that  was  derived  was  found  to  be  independent  of  the  rates. 
Since  the  rates  are  proportional  to  the  optical  power  den¬ 
sities  of  the  pulses,  the  decision  rule  will  not  be  effected 
by  signal  fluctuations. 

Ob j  ec  ti ve  of  Adding  an  AGC 

In  the  previous  section  an  optimum  decision  rule  was 
developed  for  the  detector  signal.  An  A SC  is  now  added  to  the 
optical  receiver.  The  purpose  of  the  AGC  will  be  to  amplify 
the  detector  output,  eliminate  the  signal  fluctuations  and 
minimize  the  probability  of  bit  errors.  These  objectives  will 
be  the  criteria  from  which  the  performance  of  the  AGC's  will 
be  judged. 
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III.  Receiver  Model 

In  this  chapter,  the  mathematical  models  of  the  receiver 
will  he  developed.  These  models  will  form  the  basis  for  the 
computer  simulations  and  the  analysis.  Therefore,  these 
models  must  accurately  represent  the  actual  physical  devices. 
The  model  for  the  photodetector  will  be  described  in  the  first 
section  of  this  chapter,  and  in  the  last  section,  models  for 
the  AGC  will  be  presented. 

Pho  todetec  tor  Model 

The  photodetector  model  must  accurately  describe  optical 
to  electrical  conversion  process.  To  accomplish  this  the 
photodetector  will,  at  first,  be  described  as  an  ideal 
photodetector  and  then  as  a  non-ideal  detector.  The  ideal 
photodetector  will  be  defined  as  a  detector  in  which  elec¬ 
trical  output  is  present  only  when  the  detector  is  being  il¬ 
luminated  by  an  optical  source.  In  reality,  a  small  amount  of 
electrical  output  will  be  present  even  when  no  optical  power 
is  applied;  this  will  be  considered  in  the  non-ideal  case. 

For  either  model  the  detector  will  be  followed  by  an  RC 
filter  (see  Figure  3)«  R  is  the  parallel  equivalent  of  the 
internal  resistance  losses  of  the  detector  and  an  external 
resistor,  and  C  is  the  parallel  equivalent  of  the  internal 
capacitance  losses  and  an  external  capacitor.  The  detector 


output  current  is  i(t)  and  x(t)  is  the  output  voltage  of  the 
filter. 


Figure  3*  Photodetector  Configuration 


Ideal  Pho  todetector.  An  ideal  photodetector  will  produce 
an  electrical  output  only  when  the  detector  is  being  il¬ 
luminated  by  an  optical  source.  The  photons  emitted  by  the 
source  will  strike  the  surface  of  the  detector  at  a  rate 
defined  by: 

Ap(t)  -  p ( t ) A  /  hfQ  (15) 

where 

p(t)  *  irradiance  or  incident  power  density 
A  *  effective  area  of  the  detector 
hf^  ■  energy  of  a  photon 

If  a  photon  has  sufficient  energy,  it  will  excite  an  electron 
on  the  surface  of  the  detector  into  the  conduction  band.  The 
excited  electron  will  then  be  emitted  out  of  the  detector. 


.Ik 
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This  process  will  loosely  be  defined  as  a  photon  to  electron 
conversion . 

As  more  photon  to  electron  conversions  occur,  a  detector 
output  current  can  be  described  as: 

I 

i(t)=q£5(t-t.)  (16) 

i*0  1 

where 

q  =  charge  of  an  electron 

t.  =  time  at  which  a  photon  to  electron 
conversion  occurred 

This  current  will  flow  through  the  filter  to  produce  an  output 
voltage,  x(t),  defined  by  the  convolution  integral: 

x(t)  -  /gi(T )h( t-T )  dt  (17) 

where 

_  X.  /  T>  /! 

h(t)  •  1/C  e~  =  filter  impulse  response 

or,  equivalently: 

x(t)  =*  1/C  /Qi(T)e"(t"T^/RC  dT  (18) 

If  the  period  of  observation  is  much  smaller  than  the  RC  time 
constant,  an  approximation  can  be  made  to  eqn  (18)  which 
results  in: 

x(t)  -  1/C  /Ji( t)  dt  (19) 

Substituting  in  eqn  (16)  for  the  current  and  integrating 
produces : 

I 

x(t)  -  q/C  Z  u(t-t  )  -  q/C  N(  t)  (20) 

i-0 


1  5 


r 


The  output  voltage  is  now  modelled  as  a  weighted  counting 
process,  counting  the  electrons  emitted  from  the  detector. 

Because  the  event  times  (t^'s)  of  these  photon  to  elec¬ 
tron  conversions  occur  randomly,  the  counting  process  is  Pois¬ 
son  distributed.  The  rate  of  occurrence  of  the  photon  to 
electron  conversions  is  proportional  to  the  rate  of  the 
photons  striking  the  surface  of  the  detector.  This  rate  is 
defined  as: 

X(t)  =  (np(t)A)  /  hfQ  (21) 

where 

n  *  quantum  efficiency  of  the  photon  to 
electron  conversions 

Because  the  rate  is  time  varying,  the  output  current  and 
voltage  become  inhomogeneous  Poisson  processes  in  which  the 
statistics  are  known  (Ref  8:432+)  in  terms  of  the  rate.X(t). 

Non-ideal  Pho  tod  e  tec  tor .  In  this  section  the  photodetec¬ 
tor  model  will  be  expanded  to  take  into  account  the  dark  cur¬ 
rent,  which  is  present  whether  or  not  the  detector  is  being 
illuminated.  The  dark  current  is  a  result  of  thermally  ex¬ 
cited  electrons  reaching  the  conduction  band.  Therefore,  like 
the  ideal  model,  the  dark  current  is  represented  by: 

J 

i,(  t)  *  q  Z  <5  ( t  —  t . )  (22) 

where 

t,  *  time  at  which  a  thermally  excited  electron 
J  enters  the  conduction  band 
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The  event  times  (t.'s)  also  occur  randomly  with  a  constant 

J 

rate  of  occurrence  A^  (assuming  no  significant  temperature 
variations) . 

Including  dark  current  into  the  model,  the  current  output 
becomes : 

i(t)  -  ig(t)  +  id(t)  (23) 

where 

ig(t)  =*  signal  current 
id(t)  =  dark  current 

Substituting  eqn  (16)  and  eqn  (22)  into  eqn  (23)  produces: 

I  J 

i (  t )  -  q[  I  6 ( t- 1 .  )  +  £  6 ( t- 1 .  )  ]  (24) 

i=0  1  J-0  3 

and  integrating  to  get  the  output  voltage: 

I  J 

x(t)  ■  q/C  [  Z  u(t-t  )  +  £  u(t-t  )]  (25) 

i-0  j=0  3 

Again  the  statistics  are  known  in  terms  of  the  rates  A(t)  and 

Ajj. 

The  output  voltage  given  by  eqn  (25)  will  be  the 
mathematical  model  used  to  represent  an  actual  photodetector 
output.  And  as  long  as  the  approximation  RC  is  greater  than 
the  observation  period,  this  model  will  accurately  describe 
the  photodetector. 

AOC  Model 

In  this  section  several  AGC  models  will  be  briefly  ex¬ 
amined  and  a  steady  state  response  for  a  specific  AGC  will  be 
analyzed.  For  an  analysis  of  other  AGC's  the  reader  is  refer- 


red  to  Senior  [Ref  9]*  Much  of  the  analysis  on  AGC’s  presen¬ 
ted  in  this  section  follows  the  same  approach  developed  by 
Senior. 

An  AGC  consists  of  two  major  components: 

1)  Variable  Gain  Amplifier  (VGA) 

2)  Gain  Control  Processor  (GCP) 

The  two  components  will  be  examined  separately  in  the  fol¬ 
lowing  sections.  Then  the  various  configurations  which  can  be 
derived  from  combinations  of  these  two  components  will  be  ex¬ 
amined  briefly.  In  the  last  section  a  steady  state  analysis 
will  be  conducted  for  a  specific  AGC. 


Figure  4.  VGA 


VGA.  The  VGA  is  a  three  port  amplifier  in  which  the  extra 
port  is  used  to  control  the  gain  of  the  amplifier  (see  Figure 
4).  The  VGA  is  characterized  by  its  gain  to  control  voltage 
relationship.  Virtually  every  VGA  can  be  classified  as  either 
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a  linear  VGA  where  the  gain  is  represented  as  a  function  of 


the  vc(t)  by: 

g{  t)  -  Gq  -  vc ( t)  (26) 

or  an  exponential  VGA  in  which  the  gain  is  defined  as: 


-av_( t) 

g(t)  *  Gq  e 

where 


(27) 


a  >0 

The  discussion  will  be  limited  to  these  two  VGA  models.  Other 
VGA's  are  generally  designed  to  perform  special  functions  and, 
therefore,  will  not  be  considered. 

GCP.  The  GCP  processes  some  portion  of  the  signal  to 
provide  a  control  voltage,  v^(t),  to  the  VGA.  To  keep  the 
analysis  simple  the  GCP  will  be  limited  to  the  use  of  passive 
filters,  constant  gain  amplifiers,  and  constant  biasing.  With 
this  limitation  the  output  control  voltage,  vc  (t),  easi:i-y 
defined  as: 


vCo(t)  =  A  vQi(t)  *  h ( t )  +  b  (28) 

where 

v  (t)  =*  input  voltage 
A  *  constant  gain  of  amplifier 

h(t)  ■  impulse  response  function  of  the  passive  filter 

b  *  bias  voltage 
*  *  indicates  convolution 


Basic  Configuration.  By  using  various  combinations  of 
the  two  components  (VGA  and  GCP),  numerous  AGC  models  can  be 


constructed.  The  simplest  AGC  models  are  classified  as  the 
single  loop  configurations  (see  Figure  5).  It  should  be  noted 
here  that  each  configuration  does  not  represent  one  AGC  model. 
The  configuration  is  used  to  classify  a  group  of  AGC  models; 
the  actual  AGC  model  depends  on  the  component  characteristics. 


feedforward  feedback 

Figure  5*  Single  Loop  AGC  Configurations 

The  multi-loop  configurations  contain  three  or  more  com¬ 
binations  of  the  two  AGC  components.  The  multi-loop  systems 
are  divided  into  two  sub-divisions: 

1)  Compound  Systems 

2)  Complex  Systems 

A  compound  system  is  two  or  more  single  loop  configurations 
cascaded  together.  A  complex  system  includes  everything  which 
can  not  be  classified  as  single  loop  or  compound  systems. 
Figure  6  shows  two  examples  of  complex  systems. 
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Figure  6.  Complex  AGC  Configurations 


Steady  State  Analysis .  The  steady  state  response  curve 
is  an  important  characteristic  of  an  AGC.  This  curve  can  be 
used  to  predict  an  AGC's  response  to  any  input.  Also,  for 
single  loop  conf.igura  tions ,  a  steady-state  analysis  can  be 
used  to  specify  the  AGC  parameters.  In  this  section  a  steady 
state  analysis  will  be  developed  for  a  common  AGC. 

For  the  analysis  the  VGA  will  have  an  exponential  gain 
characteristic  and  the  GCP  will  be  placed  in  a  feedback  con¬ 
figuration.  The  filter  in  the  GCP  will  have  an  impulse 
response  defined  as  h(t).  This  particular  system  is  chosen 
because  it  has  a  constant  loop  gain  which  provides  greater 
stability  (Ref  12).  An  AGC  is  considered  stable  when  the  out¬ 
put  settles  to  a  constant  value  after  an  abrupt  change  in  the 
input  (Ref  9:14). 
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Using  this  model  (see  Figure  7),  a  steady-state  analysis 
of  the  AGC  is  conducted.  The  output  of  the  AGC  is  described 
by : 

y(t)  -  x(t)  g(t)  (29) 
Substituting  in  eqn  (27)  for  the  gain  and  A  v(t)  for  the  con¬ 
trol  voltage  produces: 

y(t)  -  x(t)  Go  e-“Av(t)  (30) 
Since  the  parameters  A  and  a  have  the  same  effect  on  the  out¬ 
put,  the  parameter  a  will  be  absorbed  into  the  parameter  A. 
Thus,  the  AGC  output  is  given  by: 

y(t)  -  x ( t )  Gq  e"Av(t)  (31) 

where 

v(  t)  -  y(  t)  *  h(  t) 
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The  steady  state  response  is  obtained  by  driving  the  AGC 
with  a  constant  amplitude  input,  k^.  After  the  initial  tran¬ 
sients  decay,  a  stable  AGC  will  produce  a  constant  amplitude 
output,  kQ.  For  a  constant  output  the  filter  response  is 
defined  as: 

v(t)  -  kQ  *  h(t)  (32) 

or,  equivalently: 

v ( t )  *  kQ  *  H(0)  (33) 

where 

H(O)  -  ?{h(t)}  |f-(J 

Substituting  in  kQ  for  the  output,  k.  for  the  input,  and 
kQH(0)  for  the  control  voltage,  eqn  (31)  becomes: 

-Ak  K(0) 

kQ  -  kt  Gq  e  (34) 

Eqn  (34)  is  the  steady-state  response  of  the  AGC  and  the 
response  curve  (see  Figure  8)  is  obtained  by  plotting  k^  ver¬ 
sus  k  . 
o 

This  steady-state  response  can  serve  as  a  useful  tool  to 
obtain  the  values  of  the  AGC  parameters  Gq  and  A;  provided  the 
input-output  ranges  are  known.  For  the  analysis,  the  input  is 
assumed  to  vary  between  k^  and  k^*  The  output  is  to  be 
limited  to  the  range  kQ^  to  k^  (see  Figure  8). 


Substituting  each  pair  of  points  into  eqn  (34)  results  in: 


ki1  "  ko1  e 


AH(0)k 


ol 


/  G. 


(55) 


and 


AH(0)k 


o2 


ki2  *  ko2  6 


/  G. 


(56) 


Taking  the  ratio  of  the  above  equations  and  solving  for  A 
produces : 


a  =  m[(ki2/ki1)  (kft1/kn!>)]  /  Ch(o)  (k^-k..,)] 


ol  o2 


o2  ol 


(37) 


Once  a  value  of  A  is  determined,  0  can  be  obtained  by  using 


eqn  (35)  or  eqn  (36). 
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IV .  Computer  Simulation 

Because  of  the  non-linearities  of  the  AGC,  an  accurate 
analytical  analysis  of  the  receiver  would  be  very  difficult 
and  complex.  Or,  if  approximations  are  made  to  simplify  the 
analysis,  the  receiver  would  not  be  represented  accurately. 
Therefore,  an  empirical  analysis  was  conducted  and  the 
receiver  simulated  on  the  computer.  This  chapter  describes 

the  process  used  to  develop  a  computer  program  to  simulate  the 
receiver. 

Pho  tode  tec  tor  Algorithm 

The  photodetector  models  developed  in  the  previous  chap¬ 
ter-  easily  conform  to  computer  simulation.  The  first  section 
will  describe  the  algorithm  used  to  simulate  the  counting 
process.  The  algorithm  must  then  be  modified  to  be  able  to  be 
implemented  practically  and  still  be  able  to  accurately 
represent  the  physical  device.  The  last  section  deals  with 
the  approach  used  to  accurately  simulate  the  photodetector 
output  voltage. 

Simulation  of  Poisson  Counting  Process .  From  the 
results  of  the  previous  chapter,  the  photodetector  output 
voltage  was  mathematically  modelled  as  a  Poisson  counting 
process.  In  this  section  an  algorithm  to  simulate  the  coun¬ 
ting  process  will  be  presented. 
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i. 


For  a  Poisson  counting  process,  the  events  occur  at  ran¬ 
dom  times.  The  time  interval  between  two  successive  events 
are  also  random  and  independent  of  any  previous  time  interval. 
To  simulate  the  counting  process,  random  numbers  are  generated 
by  the  computer  to  represent  the  ranaom  event  times.  Using  a 
built-in  program,  random  numbers  which  are  uniformly 
distributed  from  0  to  1  are  generated.  The  uniform  random 
numbers  are  then  transformed  into  exponentially  distributed 
random  numbers  by  [Ref  10:62]: 

t  *  -In  (u)  /  X  (38) 

where 

t  ■  exponentially  distributed  random  variable 
u  *  uniformly  distributed  random  number  (0,1 ) 

X  *  constant  rate  for  the  occurrence  of  events 
The  random  number  generated  by  eqn  (38)  is  used  to  represent 
the  length  of  the  time  interval  between  two  successive  events. 

By  using  the  above  procedure  to  produce  the  random  time 
intervals,  a  Poisson  counting  process  is  easily  simulated. 
The  time  at  which  the  events  occur  is  defined  by: 


*i  *  *1-1  +  Ti 


(39) 


where  t^  -  arrival  time  of  the  ith  event 


i,  U 

Tj,  *  l  exponential  distributed  ;ndom  variable 


Ai  ^ ,,  , 
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POISSON  COUNTING  PROCESS 
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Figure  10.  Simulation  of  the  Poisaon  Counting  Procesi 
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(25) 


I  J 

x(t)  =*  q/C  [  i|Qu(t-ti)  +  j|Qu(t-tj)] 

Each  summation  in  eqn  (25)  represents  a  Poisson  counting 
process.  In  order  to  conduct  the  simulation,  both  Poisson 
processes  have  to  be  homogeneous.  The  second  summation  in  eqn 
(25)  is  considered  to  be  a  homogeneous  process,  but  the  first 
is  inhomogeneous  because  of  the  time-varying  optical  signal. 
Fortunately,  the  particular  optical  signal  being  studied 
enables  the  first  summation  to  be  analyzed  as  a  homogeneous 
process . 


Figure  11.  Photodetector  with  Switch  Implementation 


To  analyze  the  inhomogeneous  process  as  a  homogeneous 

process,  a  switch  is  added  to  the  photodetector  model  (see 

Figure  11).  The  clocking  period  (T  )  of  the  switch  is  half 

8 

the  bit  period  (T^)  of  the  Manchester  coded  optical  signal. 
During  the  extent  of  the  clock  period,  the  output  of  the 
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detector  will  be  observed.  If  the  switch  is  synchronized  with 
the  optical  signal,  the  power  density  of  the  received  signal 
will  remain  constant  during  a  switching  period;  and  this  will 
produce  a  homogeneous  counting  process  at  the  output.  At  the 
end  of  a  clock  period  the  switch  will  temporarily  short  cir¬ 
cuit  the  capacitor,  thus,  initializing  a  new  counting  process. 
Using  this  switching  scheme,  the  output  voltage  is  modelled  as 
a  periodic  series  of  homogeneous  Poisson  processes. 

To  implement  this  model,  the  two  counting  processes  are 
simulated  by  the  program  developed  in  the  previous  section. 
The  two  simulated  processes  are  added  and  multiplied  by  q/C  to 
produce  the  output  voltage.  The  output  is  then  sampled  at 
regular  time  increments,  At,  so  that  the  output  will  be  easily 
integrated  with  the  AGC  simulation. 

The  parameters  of  the  photodetector  simulation  are  now 
specified  to  realistically  represent  the  output  of  a  physical 
detector.  The  parameter  which  has  the  most  influence  on  an 
accurate  simulation  is  the  rate  (  X)  of  the  counting  processes. 
For  an  accurate  representation,  the  simulation  rates  are  re¬ 
quired  to  be  the  same  as  the  actual  rates  of  the  electron 
emissions.  For  an  ideal  detector,  the  emission  rates  are 
given  by: 

X  -  npA/hf o  (21) 

The  time  dependence  was  eliminated  through  the  implementation 
of  the  switch.  The  optical  power  of  sources  generally  ranges 
from  1  mw  to  10  mw  [Ref  5,  Table  l]  for  infrared  wavelengths. 


30 


Assuming  all  of  the  source's  output  power  is  detected,  the 
rate  is  calculated  to  be  an  order  of  magnitude  of  lO^/sec.  A 
rate  of  this  magnitude  is  impracticable  to  implement  directly. 

A  combination  of  two  methods  is  used  to  increase  the 
feasibility  of  the  implementation.  The  first  approach  is  a 
simple  time  scaling.  Thus,  the  rate  calculated  above  can  be 
equivalently  written  as  109/ysec  or  10^/nsec.  This  will  not 
cause  the  simulation  to  misrepresent  the  detector  because  the 
signal,  itself,  is  being  observed  in  the  scaled  time  units. 


Table  1 

Examples  of  Photodetector  Simulation  Parameters 
for  Various  Data  Rates 


Da  ta 

Rate 

Switch 

Period 

Max  Elec 
Emission 

Max 

Simula . 
Ra  te 

M 

(e'/cnt) 
_ _ > 

C 

( q/C ) M 
( V/coun  t) 

R 

(Ohms) 

1  MBit 

10  MBit 

100  MBit 

1  GBit 

10  GBit 

1  psec 

1 OOnsec 

1 Onsec 

1  nsec 

1 OOpsec 

1 09/ psec 

1 08/ 1 OOnsec 

n 

1 0  / 1 Onsec 

1 0^/ns  ec 

1  09/l OOpsec 

1000 

1000 

1000 

1000 

1000 

106 

105 

104 

105 

102 

lOpF 

2 .  OpF 

2  .OpF 

1  6nF 

50pF 

1  .  6x1 0-8 

8 . 0x1 0-9 

8.0x1 0-1 0 

1 .0x1 0-8 

3 • 2x1 0-7 

IK 

1  K 

1  OK 

1  K 

1  OK 

The  second  method  is  to  have  one  count  represent  M  electrons 
emitted  by  the  detector.  The  counting  process  is  now  weighted 
by  the  factor  M  and  the  factor  q/C.  Since  q/C  is  very  small 
in  magnitude,  the  M  can  be  easily  "absorbed"  into  the  factor 
q/c  without  a  misrepresentation  of  the  detector. 
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These  two  methods  are  applied  to  several  different  data 
rates  (see  Table  1).  The  data  rate  determines  the  clock 
period  of  the  switch.  The  electron  emission  rate  is  then  cal¬ 
culated  using  eqn  (21  )  (source  optical  power  =*  2mw)  and  scaled 
to  the  same  units  of  the  clock  period.  A  maximum  rate  for  the 
counting  process  simulation  is  determined  for  a  practical  im¬ 
plementation.  By  dividing  the  actual  electron  emission  rate 
with  the  simulation  rate,  the  factor  M  is  found.  A  capaci¬ 
tance  is  then  chosen  such  that  the  factor  (q/C)M  is  a  small 
quantity.  With  a  known  capacitance,  the  value  of  the  resistor 
is  chosen  to  produce  a  long  RC  time  constant  compared  to  the 

clock  period  (T  ).  A  long  RC  time  constant  is  required  to 
3 

support  the  approximation  that  the  detector  can  be  represented 
as  a  Poisson  counting  process.  By  using  the  two  approaches 
described  above,  a  practical  and  accurate  implementation  of 
the  photodetector  can  be  conducted. 

AGC  Algorithm 

In  this  section  an  algorithm  for  the  AGC  simulation  will 
be  presented.  The  algorithm  to  be  developed  represents  the 
single  loop  feedback  system  which  was  analyzed  in  the  previous 
chapter.  Using  this  algorithm  the  simulated  AGC  will  be 
driven  with  constant  amplitude  inputs  to  generate  a  steady- 
state  response  curve.  This  simulated  curve  is  then  compared 
with  a  theoretical  curve  to  verify  the  algorithm.  The 
validity  of  the  simulation  will  also  be  tested  using  an  am¬ 
plitude  modulated  (AN)  signal. 
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Simulation  o f  the  AGC .  In  the  preceding  chapter  an 
analysis  was  conducted  on  a  single  loop  feedback  system.  From 
the  analysis  an  equation  was  developed  to  mathematically 
describe  the  AGC.  This  mathematical  model  of  the  AGC  is: 

y  (  t )  -  x(t)  Gq  e“Av(t)  (31) 

where 

v(t)  =  y(t)  *  h(t) 

For  the  simulation  the  GCP  filter  is  specified  to  be  a  simple 
low  pass  filter.  Thus,  the  control  voltage  is  given  by  the 
convolution  integral: 

v(t)  =  1/RC  ^oy(T  )e_^  t-T  dt  (40) 

where 

H  *  filter  resistance 
C  *  filter  capacitance 

Developed  from  simple  circuit  analysis,  the  differential 
equivalent  of  eqn  (40)  is : 

C  dv  /dt  -  C y ( t )  -  v(t)]  /  R  (41) 

Both  equations  accurately  represent  the  filter,  but  eqn  (41 ) 
will  be  used  to  simulate  the  filter  because  of  its  simplicity. 

Using  numerical  methods  [Ref  4:139-190]  eqn  (41)  is  ap¬ 
proximated  by: 

[v(tkH  -  v(tk)]  /  At  -  [y(tk)  -  v(tk)3  /  RC  (42) 

or,  equivalently: 

v(tk>1)  "  y(tk)(dt/RC)  ♦  v(tk)(l*At/RC)  (43) 
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MMiifel 


where 


=  k^*1  sample  time 


At 


Sc  +  I  '  tk 


Sampling  eqn  (32)  at  the  same  rate  produces: 

-Av(tk) 

y(tfc)  -  x(tk)  G0  8 

Using  these  equations,  the  AGC  is  easily  simulated  on  the  com- 
puter.  A  block  diagram  of  the  algorithm  is  shown  in  Figure 


12. 


Figure  12.  Block  Diagram  of  the  AGC  Algorithm 


Steady-State  Response  Curve .  As  mentioned  earlier  an 
important  characteristic  of  an  AGC  is  its  steady-state 
response.  To  generate  a  response  curve,  the  simulated  AGC  is 
driven  with  constant  amplitude  inputs.  By  plotting  the  AGC 
response  (after  transients  decay)  to  these  various  inputs 
produces  a  simulated  steady-state  response  curve  (see  Figure 


3  ♦ 


13)-  This  ia  the  exact  result  which  would  be  obtained 
theoretically  for  the  same  AGC  parameters. 

Test  Signal  Response .  A  teat  signal  is  used  as  input  to 
the  AGC  to  verify  the  accuracy  of  the  simulation.  This  test 
signal  is  an  amplitude  modulated  sine  wave  defined  by: 

x(t)  *  (1  +  a  sin  2irf^t)(sin  2nfHt)  (45) 

where 

fjj  ■  carrier  frequency 

fT  *  modulation  frequency 

a  *  modulation  index 

The  AGC  will  be  designed  to  reduce  the  amplitude  variations 

caused  by  the  modulation  frequency  and  to  pass  the  carrier 

frequency  unaffected. 

Since  the  AGC  is  to  reduce  amplitude  modulations,  the  GCP 
must  track  the  envelope  of  the  signal.  The  GCP  will  consist 
of  a  full  wave  rectifier  and  a  low  pass  filter  for  envelope 

detection.  For  the  filter  to  pass  only  the  envelope  signal, 

the  RC  time  constant  is  required  to  be  within  the  range  TRef 
13:97]: 

1/fH  <<  R1C1  <<  1/fL  (46) 

Since  the  carrier  frequency  is  not  affected  by  the  AGC, 
an  approximate  analysis  can  be  conducted  by  working  with  the 
envelopes  and  neglecting  the  carrier  frequency.  The  most  ac¬ 
curate  approximate  method  [Ref  9:67]  is  to  use  the  steady- 
state  response  curve  as  shown  in  Figure  14.  The  shape  of  the 
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response  curve  is  determined  by  the  AGC  parameters  (A  and  Gq). 

The  AGC  parameters  to  produce  a  specific  response  curve 

can  be  computed  using  eqn's  (36)  and  (37).  These  equations 

require  knowledge  of  the  range  over  which  the  input  and  output 

signals  vary.  The  AGC  input  signal  (the  test  signal)  varies 

between  0.5  volts  and  1.5  volts  (for  a  modulation  index  a  = 

0.5)*  The  AGC  ouput  signal  will  be  restricted  to  the  range  - 

0.9  to  1.1  volts.  With  these  values  the  AGC  parameters  A  and 

Gq  are  calculated  to  be  4*5  and  100,  respectively. 

The  AGC  now  can  be  simulated  for  comparison  with  the 

graphical  analysis.  The  input  to  the  AGC  is  defined  by  eqn 

(45)  with  a  =  0.5,  f^  *  0.1/sec  and  f^  =  10/sec  (see  Figure 

15).  The  HC  time  constant  of  the  low  pass  filter  is  1  sec. 

The  AGC  parameters  are  A  *  6.2  and  G  *  100.  The  A  was  in- 

o 

creased  to  take  into  account  a  voltage  drop  across  the  filter. 
The  AGC  output  voltage  as  a  result  of  the  simulation  is  shown 
in  Figure  1 6 . 


Receiver  Algorithm 


The  receiver  algorithm  is  simply  the  integration  of  the 


photodetector  and  the  AGC  algor.' thms.  With  this  algorithm  the 
receiver’s  response  to  an  optical  input  can  be  simulated.  In 


this  section  a  method  will  be  developed  to  compute  the 
parameters  to  evaluate  the  performance  of  the  optical 
receiver . 


To  analyze  the  effects  of  the  AGC,  the  receiver  ia  more 


appropriately  described  by  the  first  and  second  order 
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statistics.  To  generate  the  statistics,  the  receiver  is 
simulated  several  times  with  the  same  optical  input.  After 
the  simulations,  the  mean  of  the  signal  for  a  sample  time  t^ 
is  computed  by  using  [Ref  8:245-246]: 

_  N 

a(tk)  =  [iI1ai(tk)l  /  N  (46) 

where 

s ( tfe )  3  sample  mean  of  the  signal  at  the  time  t^ 

s^(tk)  3  the  i**1  simulation  of  the  signal  at  time  t^ 

N  3  total  number  of  simulations  conducted 
The  signal  variance  is  then  obtained  by  using: 

_  N  _ 

°k2  =  a2(tk)  =  {i  =  1^a(tk)  "  3i(tk^2{  '  N  (47) 

where 

2 

3  sample  variance  of  the  3ignal  at  time  t^ 

From  these  statistics  the  effects  that  the  AGC  has  on  strength 
of  the  signal  and  the  signal  fluctuations  can  be  examined. 

To  study  the  effects  that  the  AGC  has  on  the  detection 
accuracy,  an  estimate  of  the  probability  of  a  bit  error  has  to 
be  computed.  The  decision  rule  that  minimizes  the  probability 
of  error  compares  the  voltage  (r^)  at  the  end  of  the  first 
half  of  a  bit  period  with  the  voltage  (r at  the  end  of  the 
bit  period.  This  decision  rule  can  be  applied  to  both  the 
photodetector  output  and  the  AGC  output  which  will  be  proven 
in  the  next  chapter.  A  bit  error  occurs  if  r1  <  r ^  when  a 
binary  one  is  transmitted,  or  if  r^  >  r ^  when  a  binary  zero  is 
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transmitted 


For  the  receiver  simulations,  the  probability  of  error 
can  be  computed  by  defining  a  random  variable  as  [Ref  8:265- 
266]  : 


z  . 
1 


(if  r1i 
r1i 


<  r 


2i 


>  r 


21 


1/2  If  ru.r21 
)  otherwise 


given  a  binary  one  transmitted 
given  a  binary  zero  transmitted 

(48) 


With  this  random  variable,  the  probability  of  error  is  ob¬ 
tained  by  using: 

N 

Prc(error)  =  z  *  /  N  ‘  (49) 

where 


Prc(error)  »  computed  probability  of  error 
z  =*  mean  of  the  random  variables  (z.) 

l ' 

The  variance  of  the  random  variable  is  defined  by: 
az  =  til^i-z)2]  /  » 


(50) 


By  using  this  variance  and  Tchebychef f ' s  Inequality,  the  cer¬ 
tainty  of  the  estimate  of  the  Prc(error)  can  be  quantified. 
Tchebychef  f' s  Inequality  is  given  by  [Ref  Pa  poul  i's :  259-260]  : 

Pr |  |Pr A( error)  -  Prc(error)|  <  e}  »  1  -  o ^  /  (Ne2)  (51) 

where 

PrA(error)  *  actual  probability  of  error 
The  probability  of  error  as  well  as  the  first  and  second  order 
statistics  will  be  used  to  determine  the  performance  of  the 


receiver. 


V. 


Resul ts 


In  this  section  an  implementation  of  several  optical 
receivers  will  be  conducted  and  the  effects  of  the  AGC's  ex¬ 
amined.  A  receiver  with  an  "instantaneous"  AGC  (IAGC)  will  be 
the  first  to  be  analyzed.  An  IAGC  is  a  special  classification 
of  an  AGC.  The  IAGC's  are  designed  to  be  more  sensitive  to 
input  changes,  whereas  the  AGC's  are  designed  to  respond  to 
the  slow  fluctuation  in  the  input  amplitude  (Ref  9:  81 ) .  Fol¬ 
lowing  the  IAGC  analysis  a  receiver  containing  an  AGC  will  be 
examined.  The  last  receiver  to  be  analyzed  contains  a  dual 
AGC,  which  is  (for  this  study)  an  IAGC  and  an  AGC  cascaded 
together. 

Receiver  wi th  an  I AGC 

Because  an  IAGC  is  more  sensitive  to  input  changes,  it 
should  respond  to  the  changes  caused  by  shot  noise  as  well  as 
the  signal  fluctuations.  The  effects  that  the  IAGC  has  on  the 
photodetector  output  will  be  analyzed  in  this  section.  This 
analysis  will  be  accomplished  by  computing  the  first  and 
second  order  statistics  along  with  the  probability  of  error. 

Simulation  Parameters .  To  conduct  the  analysis  the 
receiver  models  which  were  developed  in  the  previous  chapter 
will  be  used.  The  power  densities  of  the  optical  pulses  being 
received  are  arbitrarily  chosen  to  vary  over  a  20  db 
range:  0.2w/cnr  to  2.0w/cm  .  Using  eqn  (21)  the  rate  of 
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photon  to  electron  conversions  will  vary  from  10  to  10  e 
/sec  (Area»1mm^,  n-0.1,  and  f  ■  3-2x10^).  The  rate  of  ther¬ 
mally  emitted  electrons  is  generally  several  orders  magnitude 
smaller  than  the  signal  rates  and,  therefore,  was  arbitrarily 
chosen  to  be  lO^e’/sec.  With  these  counting  rates  and  with  a 
giga  bit  baud  rate,  the  simulation  parameters  for  the  detector 
are  given  in  table  1 . 

The  AGC  parameters,  A  and  Gq  can  be  determined  by  desig¬ 
ning  the  IAGC  to  limit  the  output  values  to  a  range  20mV  to 
40mV  for  inputs  that  vary  between  0.4uV  to  IO.OjjV.  The  output 
range  was  "arbitrarily”  chosen;  although  tradeoffs  were  made 
to  keep  papameters,  A  and  Gq ,  realistic.  The  input  values 
were  obtained  from  the  photodetector's  response.  The  IO.OuV 
is  the  expected  value  of  the  detector  at  the  end  of  a  clock 
period  when  the  maximum  optical  signal  is  being  received.  The 
minimum  input  value  was  determined  by  the  HAP  decision 
criteria  for  an  on-off  keying  digital  signal.  This  MAP  deci¬ 
sion  rule,  based  on  the  detector  output,  is  defined  as  (Ref 
3:216)  : 

®1 

X(Ta)  -  q/C  N(Ta)  <  q/C  >  gTs/ln( 1 g/  (52) 

where 

Ag  -  rate  of  photon  to  electron  conversions 

Aj  -  rate  of  thermally  emitted  electrons. 

Using  the  photodetector  parameters,  A^  *  10^  and  Ag  -  lO1^ 
(A  is  the  rate  for  photon  to  electron  conversions  correspon- 


ding  to  minimum  optical  power  densities)  the  optimum  decision 

threshold  is  computed  to  be  0 . 4u  V  .  Substituting  these  values 

into  egns  (36)  and  (37),  an  A  and  a  G  were  calculated  to  be 

o 

125  and  6x10^,  respectively. 

The  final  parameter  to  be  specified  is  the  RC  time  con¬ 
stant  for  the  AGC  feedback  filter.  As  long  as  the  inverse  of 

the  RC  time  constant  is  large  compared  to  the  bandwidth  of  the 
signal,  the  AGC  will  be  able  to  respond  to  rapid  changes  in 
the  input.  For  the  analysis  RC  =  O.OInsec  will  be  used. 


Analysis .  With  these  simulation  parameters,  the  effects 
of  an  IAGC  in  an  optical  receiver  can  be  studied.  For  the 
analysis  the  received  optical  signal  will  be  pre-defined  (see 
figure  17).  The  photodetector  detector  output  for  this  signal 
can  be  easily  determined  in  terms  of  the  first  and  second 
order  statistics.  Using  this  detector  output  and  the  steady- 
state  response  curve,  the  IAGC's  response  can  be  predicted 
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Figure  18.  Graphical  Analysis  of  the  TAGC  Response 
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(see  figure  18). 

The  AGC  response  generated  from  the  graphical  analysis 
indicates  that  the  IAGC  provides  control  to  the  changes  in  the 
signal.  The  relative  range  over  which  the  signal  fluctuates 
is  reduced  by  the  IAGC.  This  is  readily  apparent  by  comparing 
signal  peaks  during  the  intervals  at  which  the  optical  pulses 
were  received.  The  IAGC  will  also  reduce  the  shot  noise  ef¬ 
fects  as  evident  from  a  decrease  in  the  variance  of  the  sig¬ 
nal  . 

In  the  process  of  reducing  changes  in  the  signal,  the 
IAGC  also  distorts  the  detector  signal.  Prom  this  distorted 
signal  the  digital  information  that  was  transmitted  is  to  be 
extracted.  Before  the  addition  of  the  AGC,  the  decision  rule 
to  detect  the  information  was  defined  as 


x((2i-l)T  )  <  x((2i)T  ) 

Do 

where 

Tg  *  switch  period 

2i  T  -  iT_  =  ith  bit  period 
S  D 

x(t)  »  photodetector  output 


(55) 


Since  this  decision  rule  is  optimal  for  the  photodetector 
voltage,  applying  the  same  decision  rule  to  the  AGC  output 
would  be  appealing.  For  the  particular  AGC's  being  studied, 
the  input-output  characteristic  (steady  state  response  curve) 
is  a  mono tonically  increasing  curve  producing  a  one-to-one 
input-output  relationship.  Because  of  this  one-to-one 
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relationship,  the  same  decision  criteria  can  be  used  for  the 
AGC  output. 

The  analysis  conducted  thus  far  utilized  graphical  tech¬ 
niques  to  describe  the  effects  of  an  AGC.  Prom  the  graphic 
analysis,  an  IAGC's  response  to  the  a  photodetector  input  was 
predicted.  The  optimum  decision  rule  for  this  IAGC  response 
was  determined  to  be  the  same  criteria  that  is  used  for  the 
detector  output.  The  graphical  analysis  provides  some  insight 
as  to  what  the  IAGC  response  will  be,  but  can  not  be  used  to 
predict  the  actual  response. 

To  study  the  dynamic  effects  of  the  IAGC  requires  the  use 
of  a  computer.  Employing  the  algorithms  developed  in  the 
preceding  chapters,  a  simulation  of  the  receiver  was  conduc¬ 
ted.  The  receiver's  response  to  an  optical  input  was  obtained 
using  the  simulation  parameters  developed  in  the  first  section 
of  this  chapter.  The  statistics  of  the  photodetector's  output 
(see  figure  19)  as  a  result  of  100  simulations  accurately  ap¬ 
proximates  the  analytical  statistics.  And  the  IAGC's  response 
(see  figure  20)  to  the  detector’s  output  was  correctly  predic¬ 
ted  by  the  graphic  analysis. 

Prom  these  plots,  the  reduction  in  the  signal  variance 
may  not  be  apparent.  Therefore,  a  signal  to  noise  ratio  will 
be  computed.  The  signal  to  noise  ratio  (SNR)  is  defined  as 

SNR  -  (8k)2/s£  (54) 

where 

s^  -  sample  mean  at  time  t^ 
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Figure  19-  IAGC  Input 
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2 

3^  =  sample  variance  at  time  t^ 

And  an  improvement  factor  (IMP)  will  be  defined  aa: 

IMP  =  SNR  /SNR  (55) 

y  x 

where 

SNRy  =  SNR  of  the  AGC  output 

SNRx  *  SNR  of  the  AGC  input 

An  IMP  >  1  would  imply  the  AGC  is  reducing  the  shot  noise  ef¬ 
fects.  The  IMP  is  plotted  as  a  function  of  the  input  ensemble 

average  (see  figure  21)  which  indicates  the  IAGC  is  reducing 
the  signal  variance. 

Varying  the  AGC  parameters,  A  and  Gq,  does  not  effect  the 
general  shape  of  the  IAGC's  response.  Increasing  G^  only 
provides  greater  amplification  to  the  photodetector  output 
(see  figures  A-1  and  A -2).  Increasing  A  decreases  the  am¬ 
plification  but  provides  better  control  of  the  signal  fluctua¬ 
tions  (see  figures  A-3  and  A-4).  Varying  the  RC  time  con¬ 
stant,  on  the  other  hand,  causes  transcient  effects  to  become 

noticable  (see  figures  A-5  and  A-6).  As  the  RC  approaches  the 
switching  period,  T8>  the  IAGC  begins  to  track  the  time  or  low 
frequency  components  of  the  signal  and  becomes  less  sensitive 
to  input  changes.  When  the  reaction  time  becomes  sufficiently 
slow  (RC  ■  T  ),  the  AGC  loses  its  classification  as  an  IAGC. 

3 

As  apparent  from  the  graphical  analysis  and  through  the 
comparisions  of  the  IAGC  input-output  simulations,  the  IAGC 
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IMPROVEMENT  FACT 


Pigure  21 .  Improvement  of  SNR  by  the  IAGC 

(IMP-SNR  / SNR .  ) 

out  in 
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does  reduce  the  signal  variance.  By  reducing  the  shot  noise 
effects  the  IAGC  will  improve  the  signal  to  noise  ratio  (SNR). 
At  first  intution,  an  improved  SNR  would  indicate  that  the 
probability  of  error  would  decrease.  For  linear  systems  this 
would  be  true,  but  AGC's  have  non-linear  gain  characteristics 
which  also  shifts  the  energy  levels. 


Detector  Pr(error) 
(with  90^  certainty) 

IAGC  Parameters 

IAGC  Pr(error) 
(with  90%  certainty) 

G 

0 

■1 

RC 

0 .052+0 . 022 

2x1  o]? 
4x10 1 
6x10 1 
8x10^ 

1 0x1 05 

125 

0.01 

0.050+0.022 

0.048+0.022 

0.048+0.022 

0.047+0.022 

0.047^0.022 

6x105 

m 

0.01 

0.048+0.022 

0.048+0.022 

0.048+0.022 

0.047+0.022 

0.047^0.022 

6x1  05 

125 

0.01 

0.05 

0.1 

0.5 

1  .00 

0.048+0:022 

0.050+0.022 

0.052+0.025 

0.083+0.028 

0. 0+0.0 

Table  2.  Estimate  of  the  Probability  of  a  Bit  Error  for 
an  Initial  Signal  (binary  one)  Using  an  IAGC 
(Ag  -  5000  and  XD  -  50000) 

To  analyze  the  effects  of  the  IAGC  on  the  probability  of 
error,  again,  requires  a  simulation  of  the  receiver  on  the 
computer.  An  estimate  of  the  probability  of  a  bit  error  is 


shown  in  Table  2  for  the  first  bit  of  received  information 


(binary  one).  A  very  weak  optical  signal  and  a  large  dark 
current  were  needed  to  provide  a  non-zero  estimate  from  the 
1000  simulations.  From  the  one-to-one  input-output  relation¬ 
ship  the  graphical  analysis  predicted  that  the  IAGC  would  have 
no  effect  on  the  probability  of  error.  As  long  as  transients 
are  insignificant,  the  estimate  is  virtually  unaffected  by  the 
IAGC.  As  the  RC  time  constant  increases,  the  transients 
become  a  factor  in  the  detection  process  causing  the 
probability  of  error  to  increase. 


Figure  22.  Worst  Case  Signal 


The  influence  of  transients  on  the  detection  process  is 
more  apparent  for  the  "worst  case"  signal.  This  signal  as¬ 
sumes  a  frame  interval  of  strong  optical  pulses  are  received 
followed  by  a  frame  interval  of  pulses  from  a  weak  source  (see 
figure  22).  The  effects  of  the  IAGC  for  various  parameters  is 
shown  by  Table  3-  As  expected,  increasing  the  RC  time  con¬ 
stant  slows  the  IAGC  response,  thereby,  increasing  the 


probability  of  a  bit  error. 


Detector  Pr(error) 

IAGC  Parameters 

IAGC  Pr(error) 

(with  90^  certainty) 

G 

0 

U 

RC 

(with  90^  certainty) 

2x1  0 1 
4x1  Og 
6x10 1 
8x1  °c 

10x10’ 

1  25 

0.01 

0.062+0.024 

0.062+0.024 

0.064+0.024 

0.064+0.024 

0.064^0.024 

0.052+0.024 

75 

0.062+0.024 

6x1  05 

1  00 

1  25 

0.01 

0.063+0.024 

0.064+0.024 

1  50 

0.064+0.024 

175 

0.064^0.024 

0.01 

0.064+0.024 

6x1  0^ 

125 

0.05 

0.1 

0.067+0.025 

0.069+0.025 

0.5 

0.11 3+0.031 

1  .0 

0.204+0.040 

Table  3*  Estimate  of  Probability  of  a  Bit  Error  for 
the  "Worst  Case"  Signal  Using  an  IAGC 
(*g  ■  5000  and  =*  50000; 
for  the  previous  pulse  of  Xg=20000) 

Based  on  the  analysis  one  can  obtain  a  feel  for  the  ef¬ 
fects  of  the  IAGC  on  average  probability  of  error.  As  long  as 
the  IAGC  has  a  quick  response  to  input  charges  (RC  <<  T  )  the 
probability  of  a  bit  error  is  unaffected.  But  as  RC  in¬ 
creases,  the  bit  error  increases  causing  the  average 
probability  of  error  to  increase.  Although  the  IAGC  does  not 
improve  the  probability  of  error,  it  does  provide  excellent 
control  to  the  signal  fluctuations,  and  sufficient  amplifica¬ 
tion  of  the  signal. 
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Receiver  with  an  AGC 

The  most  common  purpose  of  the  AGC's  used  today  is  to 
control  slow-time  varying  fluctuations  that  are  super-imposed 
on  a  signal.  In  this  section  an  analysis  is  conducted  to  ex¬ 
amine  the  effects  of  a  common  AGC  in  controlling  the  optical 
signal  fluctuations. 

Simula  tion  Parameters.  For  the  analysis  the  photodetec¬ 
tor  parameters  will  he  the  same  as  those  used  in  the  previous 
section.  The  AGC  parameters  will  be  obtained  using  similiar 
methods  but  designed  to  control  the  fluctuations  and  to  pass 
the  signal  unaffected.  For  the  particular  signals  being 
analyzed,  the  fluctuations  can  be  described  as  a  change  in  the 
DC  component  or  the  time  average  of  the  signal  between  the 
frame  periods  (see  figure  23)-  Thi3  "DC"  signal  will  be 


Figure  23.  Approximation  to  Describe  Signal  Fluctuations 


used  to  compute  the  AGC  parameters  and  to  conduct  an  approx¬ 


imate  analysis  of  the  AGC  response. 


The  AGC  will  be  designed  to  limit  the  DC  output  signal  to 

the  range  to  2mV  to  4mV  for  an  input  DC  signal  that  varies 

between  0.25vV  and  25uV.  The  input  values  were  obtained  by 

taking  the  time  average  of  the  photodetector's  ensemble 

averaged  output.  The  0.25uV  corresponds  to  the  time  average 

of  the  detector’s  output  during  the  frame  period  that  a  weak 

optical  signal  was  received  and  the  25vV  corresponds  to  the 

frame  period  a  strong  optical  signal  was  received.  The  output 

range  was  arbitrarily  chosen.  Using  these  values  in  eqn  (36) 

and  eqn  (37),  the  AGC  parameters,  A  and  Gq,  are  calculated  to 
3  3 

be  2x10  and  4x10  ,  respectively.  The  final  AGC  parameter, 
the  HC  time  constant,  is  restricted  to  the  range  Tg  <  RC  <  T^.. 
This  will  insure  that  the  AGC  will  be  able  to  respond  to  the 
signal  fluctuations  and  not  respond  to  the  signal  itself. 

Analysis .  Using  the  DC  signal,  an  approximate  analysis 
can  be  conducted  through  the  graphical  techniques  (see  figure 
24).  The  result  of  the  graphical  analysis  indicates  that  the 
AGC  will  reduce  the  range  of  the  DC  variations.  Since  the  DC 
level  is  directly  proportional  to  the  average  signal;  the 
reduction  in  the  variation  of  the  DC  level  would  imply  that 
the  signal  fluctuations  will  be  reduced,  assuming  the  signal 
is  not  affected  by  the  AGC.  Because  of  the  assumptions  (the 
AGC  completely  passes  the  signal  and  the  filter  only  responds 
to  the  DC  level  of  a  frame  period),  this  graphical  analysis  is 
only  able  to  provide  a  "loose"  approximation  to  the  AGC 
response . 
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To  obtain  an  accurate  response,  the  receiver  is  implemen¬ 
ted  on  the  computer.  Using  the  optical  signal  shown  in  figure 
25,  the  photodetector's  response  (see  figure  26)  as  a  result 
of  the  simiulation  are  analogous  to  the  analytical  statistics. 
Using  an  RC  *  4.0nsec,  the  AGC  output  (see  figure  27 )  has  very 
noticable  transient  effects  which  distorts  the  photodetector 
signal . 


P(t)  (mU) 


Tt  zTf  3Tf 

Figure  25-  Input  Optical  Signal  for  a  Receiver  with  an  AGC 


The  effects  of  the  AGC  parameters  produce  similar  results 
in  the  performance  of  the  receiver  as  they  did  for  an  IAGC. 
Increasing  the  Gq  provides  greater  amplification  (see  figures 
A-7  and  A-8) ,  and  increasing  the  A  provides  better  control  of 
the  signal  fluctuations  (see  figures  A-9  and  A-10).  Like  the 
IAGC,  the  RC  time  constant  has  the  greatest  influence  on  the 
AGC  performance.  Increasing  RC  insures  the  AGC  tracks  only 
the  DC  component  of  a  frame  interval,  but  this  also  slows  the 
AGC's  response  to  the  changes  in  the  DC  level  (see  figures  A- 
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TIME  (nsec) 


Figure  27*  ACC  Output 

(G  ”4x10^,  A“2000  and  RC*4nsec) 
o 
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11  and  A-12)  introducing  transients.  Decreasing  RC  increases 
the  AGC'a  reaction  to  DC  charges  but  unfortunately  also  causes 
the  AGC  to  respond  to  the  signal.  Decreasing  RC  below  T3 
causes  the  AGC  to  be  classified  as  an  IAGC. 


De  tec  tor 

Pr( error ) 

AGC 

Pararae  ters 

AGC  Pr(error) 

(with  90? 

certainty) 

G 

0 

A 

RC 

(with  90?  certainty) 

1x10^ 

0.657+0.047 

2x10' 

0.240+0.042 

4x10' 

2000 

4.0 

0.079+0.027 

6x10,? 

0.061+0.024 

8x10’ 

0.055^0.023 

0. 0+0.0 

1  600 

0.033+0.018 

C 

1  800 

0.049+0.022 

4x1  O' 

2000 

4.0 

0.079+0.027 

2200 

0.1 58+0.036 

2400 

0.298+0.046 

'“2.0 

0.072+0.026 

4x105 

2000 

4.0 

6.0 

0.079+0.027 

0.1 56+0.036 

8.0 

0.218+0.042 

10.0 

0.241+0.043 

Table  4*  Estimate  of  the  Probability  of  a  Bit  Error  for 
the  "Worst  Case"  Signal  Using  an  AGC 
(*s  *  10000  and  -  50000; 
for  the  previous  pulse  of  Ag*20000) 

The  effects  of  the  AGC  on  probability  of  error  should 
ideally  be  neglectible,  because  an  AGC  is  to  pass  the  signal 
undistorted.  But  because  of  the  transients,  the  AGC  causes 
the  probability  of  error  to  increase  as  evident  from  Table  4 
(for  the  "worst  case"  signal  given  in  figure  22).  Note,  the 
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signal  rate  Ag  is  larger  in  this  analysis  than  in  the  IAGC 
analysis.  As  apparent  from  comparing  Tables  3  and  4,  the  AGC 
has  a  greater  detrimental  effect  on  the  average  probability  of 
error  even  though  it  is  designed  not  affect  the  signal. 

Receiver  with  £  Dua 1  AGC 

The  analysis  of  the  AGC’s  conducted  thus  far  was  limited 
simple  AGC  configurations.  In  this  section  a  simulation  of  a 
dual  AGC  (two  simple  AGC * s  cascaded  together)  will  be  conduc¬ 
ted.  The  compound  AGC  configuration  will  consist  of  an  IAGC 
followed  by  a  common  AGC.  Because  the  IAGC  has  a  quick 
response,  it  will  be  able  to  control  the  large  signal  fluctua¬ 
tions  without  affecting  the  probability  of  error.  The  AGC 
will  be  able  to  provide  further  amplification  of  the  signal 
and  provide  greater  control  of  the  fluctuations.  Since  the 
IAGC  signif i.cantly  decreases  the  dynamic  range  of  the  signal 
fluctuation,  the  transient  effect  of  the  AGC  on  the  probabili¬ 
ty  of  error  will  be  minimal. 

The  IAGC  to  be  used  in  this  compound  configuration  will 
be  the  same  as  the  IAGC  that  was  studied  in  the  first  section. 
The  IAGC  output  is  then  used  as  an  input  to  the  AGC.  The  AGC 
now  designed  to  track  the  changes  in  the  DC  component  of  the 
IAGC  output  instead  of  the  photodetector's  output.  Using  the 
same  method  as  developed  in  the  previous  section  the  AGC 
parameters  are  computed  to  be  A*0.5>  Gq**3200  and  RC*4.0nsec. 

Simulating  this  receiver  produces  an  output  shown  in 
figure  28,  where  the  detector  output  is  given  in  figure  26. 
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AGC  OUTPUT 


Figure  28.  Dual  KdC  Output 
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As  evident  from  figure  28  the  dual  AGC  is  able  to  provide  con¬ 
trol  to  the  signal  fluctuations  and  to  provide  ample  gain  to 
the  detector  output.  Since  the  dual  AGC  is  two  simple  AGC's 
cascaded  together  varying  the  parameters  should  produce  the 
same  effects  as  they  did  for  each  of  the  separate  AGC's. 
Because  of  the  time  constraints,  the  plots  of  the  dual  AGC 
ouput  for  various  parameters  could  not  be  generated. 

The  effects  that  the  dual  AGC  has  on  probability  of  error 
is  shown  in  Table  5  and  Table  6.  Through  comparisons  of  Table 
5  with  Table  4,  the  probability  of  error  is  generally  better 
than  the  AGC  alone.  This  is  expected  because  the  IAGC  is  able 
to  reduce  the  signal  fluctuations,  thereby,  reducing  detection 
errors  caused  by  the  slow  AGC  response.  Although  comparing 
Table  6  with  Table  3,  the  AGC's  slow  response  does  still  have 
a  detrimental  effect  on  detection  accuracy. 
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IAGC  Parameters 


(x  -10000  and 

xJ-50000; 

for  previous 
pulse  of 
Xg-20000) 


AGC  Parameters  1  Dual  AGC  Pr(error) 
RC  I  with  90i  certainty 


0.0270.0162 

0.008+0.0089 

0.0030.0055 

0.0030.0055 

0.003+0.0055 


0.0030.0055 

0.0030.0055 

0.0030.0055 

0.0070.0083 

0.0130.0113 


0.0030.0055 

0.0060.0077 

0.0080.0089 

0.1630.0369 

0.5500.0497 


0.0270.0122 

0.0070.0083 

0.0030.0055 

0.0030.0055 

0.0030.0055 


0.0000.0000 
0.0030.0055 
0. 0030.0055 
0.0140.0117 
0.0570.0232 


0.0430.0203 
0.0030.0055 
0.0000.0000 
0.0000.0000 
0.000 +0.0000 
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Detector  Pr(error) 
with  90%  certainty 

IAGC  Parameters 

AGC  Parameters 

Dual  AGC  Pr(error) 
with  certainty 

Go 

A 

RC 

G0 

■ 

RC 

0.063+0.024 

2x1  o’ 
4x10 1 
6x10  \ 
8x1  o’ 
10x1 O’ 

125 

0.01 

3200 

0.5 

1 

0.309+0.0462 

0.169+0.0375 

0.103+0.0304 

0.070+0.0255 

0.055  <>.0228 

6x1 05 

75 

100 

125 

150 

175 

0.01 

3200 

1 

1 

0.053+0.0224 

0.065+0.0247 

0.103+0.0304 

0.155+0.0362 

0. 213*0.0409 

6x1  o’ 

125 

RHgn 

jjjsl 

Rllgi 

3200 

0.5 

0.103+0.0304 

0.107+0.0309 

0.135+0.0342 

0.396+0.0489 

0.782+0.0413 

6x1 05 

125 

0.01 

2400 

2800 

3200 

3600 

4000 

0.5 

1 

0.245+0.0430 

0.158+0.0346 

0.103+0.0304 

0.069+0.0253 

0.055*0.0228 

(1,5=5000  and 

A “=50000; 
for  previous 
pulse  of 
lg=20000) 

6x1 05 

125 

0.01 

3200 

0.40 

0.45 

0.50 

0.55 

0.60 

1 

0.013+0.0113 

0.040+0.0196 

0.103+0.0346 

0.229+0.0420 

0.456+0.0498 

6x1 05 

125 

0.01 

3200 

1 

9 

0.330+0.0471 

0.103+0.0304 

0.039+0.0194 

0.021+0.0143 

0.017+0.0129 

Table  6.  Estimate  of  the  Probability  of  a  Bit  Error  for 
the  "Worst  Case"  Signal  Using  a  Dual  AGC 


VI.  Conclusions  and  Recommendations 


Conclusions . 

The  purpose  of  this  paper  was  to  analyze  the  effects  of 
an  AGC  in  an  optical  receiver.  The  particular  optical  signal, 
which  was  studied,  was  a  Manchester  coded  digital  signal  with 
varying  optical  power  densities.  The  optical  signal  was  con¬ 
verted  into  an  electrical  signal  by  a  photodetector.  Inserted 
after  this  detector  was  an  AGC  with  the  intent  of  improving 
the  performance  of  the  receiver.  The  criteria  to  assess  the 
AGC's  effects  on  the  receiver's  performance  were: 

1)  amplification  of  the  detector  signal 

2)  control  of  the  signal  fluctuations 

3)  improve  the  detection  accuracy. 

Using  this  criteria  three  specific  AGC's  were  analyzed. 

From  the  analysis,  an  AGC  which  has  a  quick  response  was 
determined  to  be  the  most  suitable  for  the  optical  receiver. 
Although  a  quick  AGC  response  causes  distortion  of  the 
photodetector  output,  the  decision  criteria  was  not  effected 
by  the  addition  of  the  AGC's.  Of  the  three  AGC’s  examined, 
the  IAGC  had  the  fastest  reaction  to  input  changes  (as  long  as 
RC<<T  ).  The  IAGC  was  able  to  provide  amplification  of  the 
signal  and  to  control  signal  fluctuations,  but  it  had  no  sig¬ 
nificant  effect  on  the  probability  of  error.  Because  of  their 
slow  response,  the  other  AGC's  that  were  studied  caused  an 
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increase  in  probability  of  error.  The  IAGC  was  able  to 
provide  amplification  of  the  signal  and  control  signal  fluc¬ 
tuations,  but  it  had  no  significant  effect  on  the  probability 
of  error.  Thus,  an  AGC  in  an  optimal  receiver  is  not  capable 
of  improving  the  probability  of  error.  Although  this  is  not 
the  result  which  is  desired,  it  is  an  intuitively  pleasing 
result.  Having  a  device  (which  controls  the  signal  amplitude) 
improve  the  accuracy  of  PPM  detection  would  not  have  been  an¬ 
ticipated. 

Although  the  AGC  will  not  improve  the  detection  accuracy, 
it  can  be  used  to  amplify  the  signal  and  to  control  the  signal 
fluctuations.  For  signals  with  a  wide  dynamic  range  of  fluc¬ 
tuations,  using  an  AGC  in  the  receiver  may  be  more  applicable 
than  a  constant  gain  amplifier.  Since  the  AGC  had  non-linear 
gain  characteristics,  it  can  be  use  to  provide  ample  gain  to  a 
weak  signal  to  simplify  the  detection  process.  And  during  the 
reception  of  a  strong  signal  the  AGC  will  automatically  reduce 
the  gain  to  prevent  saturation  of  the  decision  circuitry. 
Therefore,  circuit  limitations  may  require  the  use  of  an  AGC, 
even  though  an  AGC  will  not  improve  the  probability  of  error. 

Recommendations 

Because  of  the  numerous  possible  AGC  configurations  and 
time  constraints,  only  three  AGC's  were  examined,  therefore 
the  analysis  could  be  expanded  to  investigate  other  AGC 
models.  Since  the  AGC's  can  not  improve  the  detection  ac¬ 
curacy,  the  study  should  be  focused  on  minimizing  an  increase 


in  the  probability  of  error,  while  maintaining  that  the  AGC 
still  provide  sufficient  amplification  of  the  signal  and  con¬ 
trol  of  the  fluctuations. 

Other  photodetector  models  could  also  be  analyzed.  In 
this  paper  the  photodetector  was  modelled  as  a  weighted  Pois¬ 
son  counting  process.  This  required  a  large  RC  time  constant 
(compared  to  the  signal  bandwidth)  and  the  implementation  of  a 
switch.  Using  a  large  RC  time  constant  caused  the  signal  in¬ 
formation  to  be  lost  and  an  impracticable  switch  had  to  be 
implemented  to  obtain  the  information  sent.  Therefore,  a  more 
accurate  photodetector  model  could  be  investigated  by 
eliminating  the  switch,  which  would  require  the  RC  time  con¬ 
stant  in  the  detector  model  to  be  reduced. 

The  analysis  could  al30  be  expanded  to  examine  other  op¬ 
tical  signals  and  channels.  A  study  of  optical  signals  which 
are  transmitted  through  the  atmosphere  may  be  more  adaptable 
to  AGC  control.  The  atmospheric  turbulance  will  cause  the 
signal  to  incur  continuous  time  varying  fluctuations,  unlike 
the  discrete  fluctuations  caused  in  the  fiber  optic  system 
which  was  analyzed.  With  the  continuous  time-varying  fluctua¬ 
tions  the  AGC  will  not  be  required  to  have  a  quick  response. 
This  will  enable  the  use  of  AGC's  which  are  commonly  used 
today . 


70 


in  the  probability  of  error,  while  maintaining  that  the  AGC 
still  provide  sufficient  amplification  of  the  signal  and  con¬ 
trol  of  the  fluctuations. 
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switch.  Using  a  large  RC  time  constant  caused  the  signal  in¬ 
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accurate  photodetector  model  could  be  investigated  by 
eliminating  the  switch,  which  would  require  the  RC  time  con¬ 
stant  in  the  detector  model  to  be  reduced. 

The  analysis  could  also  be  expanded  to  examine  other  op¬ 
tical  signals  and  channels.  A  study  of  optical  signals  which 
are  transmitted  through  the  atmosphere  may  be  more  adaptable 
to  AGC  control.  The  atmospheric  turbulance  will  cause  the 
signal  to  incur  continuous  time  varying  fluctuations,  unlike 
the  discrete  fluctuations  caused  in  the  fiber  optic  system 
which  was  analyzed.  With  the  continuous  time-varying  fluctua¬ 
tions  the  AGC  will  not  be  required  to  have  a  quick  response. 
This  will  enable  the  use  of  AGC's  which  are  commonly  used 
today . 
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Appendix  A 

Graphs  of  Various  AGC  Responses 

This  appendix  contains  plots  of  the  AGC  response  to  a 
given  optical  input  for  various  AGC  parameters.  The  plots  are 
presented  in  pairs  in  which  one  AGC  parameter  is  varied 
between  the  pair.  This  will  enable  the  reader  to  compare  the 
effects  that  the  parameter  has  on  the  photodetector  output. 

The  photodetector  output  for  the  first  six  graphs  is 
given  in  Figure  19  and  the  last  six  in  figure  26.  Each  graph 
contains  a  plot  of  the  mean  +_  100  standard  deviations,  which 
is  a  result  of  100  simulations. 
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AGC  OUTPUT  VS  TIME 
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Figure  A- 3.  IAGC  with  A«75 


(G  *6x10  and  RC*0.01nsec) 
o 


75 


AGC  OUTPUT  VS  TIME 
(RECEIVER  OUTPUT) 


Figure  A-4.  IAGC  Output  with  A*150 


(Go»6x107  and  RC*0 . 0 1  ns ec ) 
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AGC  OUTPUT  VS  TIME 


Figure  A-4.  IAGC  Output  with  A-150 
(Goa6x10^  and  RC*0.01nsec) 
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Figure  A-5«  IAGC  Output  with  RC»0. Insec 
(0o*6x10^  and  A-125) 
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Figure  A-6.  IAGC  Output  with  RC*1 .Onaec 
(0o-6x105  and  A-125) 
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AGC  OUTPUT  VS  TIME 
(RECEIVER  OUTPUT) 

(mean  +  100  standard  deviations) 


Figure  A-7.  AGC  Output  wit!'.  Gq»1x105 
(A*2000  and  HC*4.0nssc) 
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AGC  OUTPUT  VS  TIME 
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Figure  A-8.  AGC  Output  with  Ga*8x10 


(A-2000  and  RC*4.0nsec) 
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Figure  A-9.  ACC  Output  with  A*1800 

(C  and  RC-4.0naec) 
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Appendix  B 


Computer  Program 

This  appendix  contains  the  program  which  was  used  to 
generate  the  data  for  the  plots  and  the  tables  provided  in 
this  report.  The  program  is  written  in  PORTRAlf  IV  programming 
language  and  implemented  on  the  CSC  6600  computer  system.  The 
uniform  random  number  generator  subroutine  is  not  presented 
here.  This  routine  was  obtained  from  the  INSL  library  which 
was  available  at  AFIT. 
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100-  PROGRAM  DUAL ( INPUT , OUTPUT , TAPE5- INPUT , T APE6- INPUT , T APE1 , 

1 1 0-  *T APE2 , TAPE3 , TAPE4 , TAPE7 , TAPE8) 

120-  DIMENSION  Xl(lOO),Yl(lOO),Vl(lOO),Tl(lOO),T2(lOO),R(2),XLAMDA(3l) 

130-  *,YT(lOO),XT{lOO),Y2(lOO),V2(lOO) 

140-  DOUBLE  PRECISION  DSEED 

150-  INTEGER  C0UNT1 .C0UNT2 

160-C 

170-C 

180-C 

1 90-C  ***  INITIAL  SEED  FOR  RANDOM  NUMBER  GENERATION  *** 

200-C 

210-  DSEED-92751 

220-C 

230-C 

240-C  ***  NUMBER  OP  SIMULATIONS  OP  A  GIVEN  OPTICAL  SIGNAL  *** 

250-C 

260-  NRUN-100 

270-C 

280-C 

2 90-C  ***  INPUT  OPTICAL  SIGNAL  GIVEN  IN  TERMS  OP  THE  SIMULATION  RATES  ** 

300-C  ***  THESE  RATES  WERE  ARBITRARILY  CHOSEN  BY  THE  PROGRAMER  *** 

31 0-C 

320-  XLAMDA(1)=0.0 

330-  XLAMDA(2)-100.0 

340-  XLAMDA(3)-100.0 

350-  XLAMDA(4)-0.0 

360-  XLAMDA(5)-100.0 

370-  XLAMDA(6)-0.0 

380-  XLAMDA(7) -100.0 

390-  XLAMDA(8)-0.0 

400-  XLAMDA(9)-0.0 

410-  XLAMDA(10)-100.0 

420-  XLAMDA(ll)-800.0 

430-  XLAMDA(l2)-0.0 

440-  XLAMDA( 13) -800.0 

450-  XLAMDA(l4)-0.0 

460-  XLAMDA( 1 5 ) -0 . 0 

470-  XLAMDA( 16) -800.0 

480-  XLAMDA(l7)-0.0 

490-  XLAMDA( 18) -800.0 

500-  XLAMDA(l9)-0.0 

510-  XLAMDA( 20) -800.0 

520-  XLAMDA( 21 ) -200 . 0 

530-  XLAMDA(22)-0.0 

540-  XLAMDA(23)-0.0 

550-  XUMDA(  24) -200.0 

560-  XLAMDA(25)-0.0 

570-  XLAHDA( 26) -200.0 

580-  XLAMDA(27)-0.0 

590-  XLAMDA(28)-200.0 

600-  XLAMDA(29)-0.0 

610-  XLAMDA(30)-200.0 


36 


«* 


630-C 
640-C 
650- 
660-C 
670-C 
680-C 
690-C 
700-C 
710-C 
720- 
730-C 
740-C 
750-C 
760-C 
770-C 
780-C 
790-C 
800-C 
81 0-C 
820-C 
830- 
840- 
850- 
860-C 
870-C 
880-C 
890-C 
900- 
910- 
920- 
930-C 
940-C 
950-C 
960-C 
970-C 
980-C 
990-C 
1000-C 
101 0-C 
1020-C 
1030-C 
1040-C 
1050-C 
1060- 
1070-C 
1080-C 
1090-C 
1100- 
111 0-C 
1120-C 
1130-C 
1140-C 
1150- 


***  SIMULATION  RATE  OF  THE  THERMALLY  EMITTED  ELECTRONS  *** 
DARK-10.0 


***  (Q/C)*M  -  C0NST1  -  WEIGHTING  FACTOR  FOR  THE  COUNTING  PROCESS 
TO  PRODUCE  AN  ACCURATE  REPRESENTATION  OF  THE  DETECTOR  *** 

C0NST1 -1 .OE-8 


***  AGC  PARAMETERS  *** 

MM HM»  M  M  M  M 

HR  Vi RR WWwWWwWWwWwWwwWw 


***  AGC  PARAMETERS  FOR  THE  FIRST  AGC  (IAGC)  *** 

RC-0.01 

G-600000.0 

A-125-0 


***  AGC  PARAMETERS  FOR  THE  SECOND  AGC  (COMMON  AGC)  *** 

RC2-10.0 

A2-0.5 

G2-3200.0 


TIME  PARAMETERS  (UNITS  ARE  RELATED  TO  THE  BAUD  RATE) 


***  PERIOD  OF  THE  SWITCH  (CLOCKING  PERIOD)  *** 
***  THIS  PERIOD  -  1/2  THE  BIT  PERIOD  *** 

PERIOD-1 .0 

***  COMPUTER  SAMPLING  INTERVAL  *** 

DELTA-0.001 


***  NUMBER  OF  CLOCK  PERIODS  TO  SIMULATED  *** 
HPER-4 


•  • 
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***  NUMBER  OP  SAMPLES  PER  CLOCK  PERIOD  *** 
NSAMP-  INT  ( ( PER  IOD-t-O .  OOOOOOI )  /DELTA  ) 


1 1 70-C 
1180*0 
1190-C 
1200- 
1210-C 
1220-C 

1230-C  ***  THE  J-TH  SAMPLE  AT  WHICH  THE  STASTICS  ARE  COMPUTED  »** 

1240-C 

1250-  JSAMP-50 

1 260-C 

1270-C 

1280-C  ***  THE  BIT  PERIOD  AT  VTHICH  PROB(ERROR)  IS  COMPUTED  *** 

1290-C 

1300-C  CLOCK  PERIOD  FOR  THE  FIRST  HALF  OF  THE  BIT  PERIOD 
1310-  NPER1-5 

1320-C  CLOCK  PERIOD  FOR  THE  SECOND  HALF  OF  THE  BIT  PERIOD 

1330-  NPER2-6 

1340-C 

1350-C 

1360-C 

1370-C 

1380-C 

]  390=*C  <HHHHW<<HHHHHHHHHHHHHHHHHHHK< 

1400-C  ***  INITIALIZE  PARAMETERS  *** 

1410-C  <<<<<<<<<<<<<<<<<<<<<<<<<<<  <  < 

1420-C 

1430-C 

1440-  DO  10  K-1  ,NRUN 

1450-  T1 (K)-O.O 

1460-  T2(K)-0.0 

1470-  XI (K)-O.O 

1480-  V2(K)-0.0 

1490-  VI (K)-O.O 

1500-10  CONTINUE 

151 O—C 
1520-C 
1 530-C 
1 540-C 

1550-C  ***  BEGIN  SIMULATION  *** 

1560-C 

1570-  DO  400  1-1 ,NPER 

1 580-C 

1590-C 

1600-  JTH-1 

1610-  COUNT 1 -0 

1620-  C0UNT2-O 

1630-C 
1640-C 

1650-C  ***  SIMULATE  ONE  CLOCK  PERIOD  *** 

1660-C 

1670-  DO  3 00  J-1.NSAM P 

1680-C 

1 690-  TIME-DELTA*FLOAT( J) 


1710-C 
1720-C 
1730-C 
1740- 
1760-C 
1770-C 
1780-C 
1 790-C 
1800-C 
1810- 
1820- 
1830-C 
1840-C 
1850-C 
1860-C 
1870-20 
1880-C 
1890- 
1900-C 
1910- 
1920-C 
1930-C 
1940- 
1 950-C 
1 960-C 
1970- 
1980-C 
1990- 
2000- 
2010-30 
2020-C 
2030- 
2040-C 
2050-C 
2060-C 
2070-C 
2080-C 
2090-C 
2100-C 
2110-C 
2120-C 
2130- 
2140- 
21 50-40 
2160- 
2170- 
2180- 
2190- 
2200- 
2210- 
2220-50 
2230- 


***  COMPUTE  THE  DETECTOR  VOLTAGE,  AGC  OUTPUT  VOLTAGE,  AGC 
CONTROL  VOLTAGE  FOR  EACH  RUN  *** 

DO  100  K-1.NRUN 


***  GENERATION  OF  A  COUNTING  PROCESS  TO  SIMULATE  AN  IDEAL  DETECTOR 

IF  (XLAMDA(I).EQ.O.O)  GOTO  30 
COUNT 1-0 

T1(K)  IS  THE  TIME  AT  WHICH  M  ELECTRONS  ARE  EMITTED  FROM  THE 
IDEAL  PHOTODECTOR  AS  A  RESULT  OF  THE  INCIDENT  OPTICAL  SIGNAL 

if(ti (k) .gt.time)  GOTO  30 

GGUBS  -  IMSL  UNIFORM  NUMBER  GENERATOR  SUBROUTINE 
CALL  GGUBS ( DSEE D , 2 , R ) 

GENERATE  A  NEW  SEED 
DSEED»R( 2)*1 23457.0 

CONVERT  UNIFORM  NUMBER  INTO  A  EXPONENTIAL  RANDOM  NUMBER 
TAU1 --ALOG(R( 1 ) ) /XLAMDA( I ) 

INCREMENT  COUNT;  ONE  COUNT  IS  EQUIVALENT  TO  M  PHOTON  TO 
ELECTRON  CONVERSIONS 
C0UNT1 -C0UNT1 *1 

TIME  AT  WHICH  M  PHOTON  TO  ELECTRON  CONVERSIONS  OCCUR 
T1 (K)-T1 (K)+TAU1 
GOTO  20 
CONTINUE 

INITIALLY  ONE  COUNT  TO  MANY  IS  GENERATED 
if((j.eq.i ) .and. (counti .ne.o})  C0UNT1 -C0UNT1  -1 


GENERATION  OF  COUNTING  PROCESS  TO  ADD  DARK  CURRENT  EFFECTS  *** 


THESE  PARAMETERS  ARE  SIMILAR  TO  THOSE  DESCRIBED  IN  THE 
IDEAL  SIMULATION 

IP(DARK.EQ.O.O)  GOTO  50 
C0UNT2-0 

IF (T2(K). GT.TIME)  GOTO  50 
CALL  GGUBS(DSEED,2,R) 

DSEED-R(2)*123457.0 
TAU2--AL0G( R( 1 ) ) /DARK 
C0UNT2-C0UNT2+1 
T2 ( K) -T2 ( K) +TAU2 
GO  TO  40 
CONTINUE 

IF(( J.EQ.1). AND. (C0UNT2. NE.O))  C0UNT2-C0UNT2-1 
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2240-C 
2250-C 
2260-C 
2270-C 
2280- 
2290-C 
2300-C 
231 0-C 
2320-C 
2330-C 
2340-C 
2350-C 
2360-C 
237 0-C 
2380- 
2390-C 
2400- 
241 0-C 
2420- 
2430-C 
2440- 
2450-100 
2460-C 
2480-C 
2490-C 
2500-C 
251 0-C 
2520-C 
2530-C 
2540-C 
2550-C 
2560-C 
2570-C 
2580- 
2590- 
2600- 
2610- 
2620- 
2630- 
2640- 
2650- 
2660-C 
2670-C 
2690-C 
2700-C 
2710-C 
2720- 
2730- 
2740- 
2750-60 
2760-C 
2770- 
27 80- 


***  PHOTODETECTOR  SIMULATION  *** 

MHHHHHHHHHHHHHHHt  WHHHHHHHHHHHUK 


XI  (K) -C0NST1 *PL0 AT ( COUNT 1 +C0UNT2 )+X1 (K) 


***  DUAL  ACC  SIMULATION  *** 


I ACC  OUTPUT  VOLTAGE 

Y1 (K)=X1 (K)*G*EXP(-(A*V1 (K))) 

IAGC  CONTROL  VOLTAGE 

VI (K)-(Y1 (K)-V1 (K))*DELTA/RC+V1 (K) 

AGC  OUTPUT  VOLTAGE 

Y2(K)«Y1 (K)*G2*EXP(-(A2*V2(K))) 

AGC  CONTROL  VOLTAGE 

V2(K)-(Y2(K)-V2(K))*DELTA/RC2+V2(K) 

CONTINUE 


WTHrW  NVH  WWWwWWIrWwWTNII  ■  II  ■  WWW 


***  COMPUTE  STATISTICS  *** 


CALCULATE  STATISTIS  FOR  EVERY  J-TH  SAMPLE 

IF ( JTH . NE . JSAMP )  GOTO  200 

JTH-0 

SUM1-0.0 

SUM2-0.0 

SUM3-0.0 

SUM4-0.0 

SUM5-0.0 

SUM6-0.0 


***  COMPUTE  MEAN  *** 

DO  60  K-I.NRUN 
SUM1 -SUM1 +X1 (K) 
SUM2-SUK2+Y2(K) 
CONTINUE 

XMN-SUM1 /FLOAT(NRUN) 
YMN-SUM2/FL0AT( NRUN) 
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2800-C 
281 0-C 
2820- 
2830- 
2840- 
2850-70 
2860-C 
2870- 
2880- 
2890-C 
2900-C 
291 0-C 
2920-C 
293 0-C 
2940-C 
2950- 
2960- 
2970-C 
2980-C 
2990-C 
3000-C 
301 0-C 
3020- 
3030-C 
3040-C 
3050-C 
306 0-C 
3070-C 
3080- C 
3090-C 
31  CO- 
31 10- 
3120-C 
3130-C 
3140-200 
3150-C 
3160- 
3170-C 
3180-C 
3190-300 
3230-C 
3240-C 
3250-C 
3260-C 
3270-C 
3280-C 
3290- 
3300- 
3310- 
3320- 
3330-71 
3340-C 
3350-72 


***  COMPUTE  VARIANCE  **# 

DO  70  K-1.NRUN 
SUM3-SUM3+(XMN-X1 (K))**2 
SUM4-SUM4 + ( YMN -Y 1 (K) )**2 
CONTINUE 

XVAR-SUM3/FL0AT(NRUN) 
YVAR-SUM4/FL0AT ( NRUN ) 


***  COMPUTE  SIGNAL  TO  NIOSE  RATIOS  *** 

SNRX-XMN**2/XVAR 

SNRY-YMN**2/YVAR 

***  IMPROVEMENT  FACTOR  *** 
IMP-SNRY/SNRX 


***  WRITE  ONTO  TAPE  NEEDED  DATA  *** 

WRITE( 1 ,*)  XMN, XVAR 
WRITE(7,*)  YMN.YVAR 

CONTINUE 

IP( JSAMP.GT. 1 )  JTH-JTH+1 
CONTINUE 

***  COMPUTE  THE  PROBABILITY  OF  ERROR  *** 


IF(NPER.NE.NPERI)  GOTO  72 
DO  71  JK-1.NRUN 
XT(JK)-X1(JK) 
YT(JK)-Y2(JK) 

CONTINUE 

CONTINUE 
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3390-C  ***  THIS  ROUTINE  WILL  ONLY  COMPUTE  THE  PR(ERROR) 

3400-C  FOR  A  BINARY  "ONE"  *** 

341 O-C 

3420-  IF( NPER . NE . NPER2 )  GOTO  75 

3430-  SUM1-0.0 

3440-  SUM2-0.0 

3450-  SUM3-0.0 

3460-  SUM4-0.0 

3470-C 

3480-C 

3490-  DO  74  JK-1 ,NRUN 

3500-  IF(X1(JK).GT.XT(JK))  SUM1-SUM1 +1 .0 

3510-  IF(X1(JK).EQ.XT(JK))  SUM2-SUM2+0.5 

3520-  IF(Y2(JK).GT.YT(JK))  SUM3-SUM3+1 .0 

3530-  IF(Y2(JK).EQ.YT(JK))  SUM4-SUM4+0.5 

3540-74  CONTINUE 

3560-C 

3570-C 

3580-C  COMPUTE  A  PROBABILITY  OF  ERROR  FOR  THE  INPUT 

3590-  PRX-(SUM1+SUM2)/FL0AT(NRUN) 

3600-C 

3610-C  COMPUTE  A  PROBABILITY  OF  ERROR  FOR  THE  OUTPUT 

3620-  PRY=(SUM3+SUM4)/FL0AT(NRUN) 

3630-C 

3650-C 

3660-C 

3670-C  COMPUTE  THE  VARIANCE  OF  THE  INPUT  PR( ERROR) 

3680-  PRXVAR-(SUM1*(l . O-PRX ) **2  +SUM2* ( 0 . 5-PRX ) **2 ) /FLOAT ( NRUN) 

3690-  PRXVAR-PRXVAR+ ( 1 . 0-PRX)*PRX**2 

3700-C 

3710-C  COMPUTE  THE  VARIANCE  OF  THE  OUTPUT  PR( ERROR) 

3720-  PRYV  AR- ( SUM3* ( 1 . O-PRY ) **2+SUM4# (0.5- PRY ) **2 ) /FLO AT ( NR  UN) 

3730-  PRYVAR-PRYVAR+(l .0-PRY)*PRY**2 

3740-C 

3750-C 

3760-75  CONTINUE 

3770-C 

3790-C 

380O-C 

3810-C  ***  RE-INITIALIZE  PARAMETERS  *** 

3820-C 

3830-  DO  80  K-1 ,NRUN 

3840-  XI (K) -0.0 

3850-  T1 (K)-O.O 

3860-  T2(K)-0.0 

3870-80  CONTINUE 

3880-C 
3890-C 

3900-400  CONTINUE 

391 0-C 

3950-  STOP 

3960-  END 

3970-*B0R 
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This  paper  is  to  analyse  the  effects  of  an  automatic  gain 
control  (ACC)  in  an  optical  receiver  system.  The  intent  of 

'using  an  AGC  is  to  improve  the  performance  of  the  receiver. 

The  received  signal  is  Manchester  encoded  and  contains 
fluctuations  caused  by  transmission  attenuation.  In  the  con¬ 

version  from  an  optical  signal  to  an  electrical  signal,  the 
photodetector  produces  shot  noise.  The  purpose  of  the  AGC 
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amplify  the  photodetector  signal,  control  signal 
,  reduce  shot  noise,  and  improve  detection  ac- 

this  criteria,  +  v>s  effects  of  the  AGC  in  the 
examined.  But  because  of  the  non-linearities 
ith  an  AGC,  an  analytical  analysis  can  not  be  con- 
refore,  the  receiver  was  simulated  on  the  corn- 
paper  describes  the  models  and  algorithms  used  to 
receiver.  Through  the  simulations,  the  receiver's 
with  an  AGC  is  compared  to  the  performance  without 

of  the  limited  time  for  the  analysis,  only  three 
*s  were  examined.  The  first  AGC  studied  is  clas- 
instantaneous  AGC  (IAGC)  followed  by  an  analysis 
AGC.  The  final  AGC  to  be  investigated  is  a  dual 
an  AGC  cascaded  to  an  IAGC. 
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